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Abstract 

The present thesis aims at an extension of the canonical formalism of Arnowitt, Deser, and Misner 
from self-gravitating point-masses to objects with spin. This would allow interesting applications, e.g., 
within the post-Newtonian (PN) approximation. The extension succeeded via an action approach to 
linear order in the single spins of the objects without restriction to any further approximation. An order- 
by-order construction within the PN approximation is possible and performed to the formal 3.5PN order 
as a verification. In principle both approaches are applicable to higher orders in spin. The PN next-to- 
leading order spin(l)-spin(l) level was tackled, modeling the spin-induced quadrupole deformation by a 
single parameter. All spin-dependent Hamiltonians for rapidly rotating bodies up to and including 3PN 
are calculated. 
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1 Introduction 

Though general relativity has seen and passed many experimental tests, one of its most fascinating predic- 
tions, namely gravitational waves, has not been observed directly. However, observations of certain binary 
pulsar signals are in good agreement with the energy loss predicted by general relativity due to gravitational 
waves, see, e.g., [1]. This indirect observation of gravitational waves originates from Hulse and Taylor (first 
found for the binary pulsar PSR B1913+16) and was awarded the Nobel Prize in 1993. Nowadays there 
is less doubt that gravitational waves exist, and one aims at a direct observation with assiduous efforts, 
both by experiments on Earth, e.g., LIGO, VIRGO, GEO 600, and by the future space mission LISA [2]. 
The direct measurement of gravitational waves is not only interesting, but would furthermore open up an 
entirely new spectrum for astronomical observations. Such gravitational wave astronomy is expected to 
have great impact on astrophysics and fundamental physics [2], possibly starting a new era in these fields. 

Beside the experimental challenge of measuring extraordinarily small relative changes in length (< 
10 -21 detectable by now) there are important problems to be solved on the theoretical side in order to 
successfully establish the new field of gravitational wave astronomy. The theoretical challenge lies within 
the area of data analysis, for both the noise dominated [3] and signal dominated [4] cases. An accurate un- 
derstanding and knowledge of the expected gravitational wave signals is a key ingredient to allow faithful 
astronomical or astrophysical statements from the data analysis process. An appealing source for gravi- 
tational waves is the inspiral and merger of two compact objects, like black holes and neutron stars. The 
advantage of this kind of source is its quite periodic behavior, which can be studied over long periods of 
time. However, minute changes in frequency and amplitude of the gravitational waves need to be predicted 
in an accurate way. While fully numerical methods are ideal to study the very late inspiral (or plunge) 
and merger phases of compact objects, the post-Newtonian approximation to general relativity provides 
a good analytic handle on the inspiral phase and can give accurate predictions over many orbits. The 
post-Newtonian approximation was pushed to high orders for nonspinning objects, see, e.g., [5], and it is 
desirable to catch up to these orders for the spinning case. 

A successful and efficient way to calculate the conservative part of the dynamics of two compact objects 
within the post-Newtonian approximation is based on the canonical formalism of Arnowitt, Deser, and 
Misner (ADM). However, this formalism has been coupled so far to nonspinning point-like objects only. 
The main goal of the present thesis is to extend this coupling to spinning objects. Not only this is useful 
for subsequent applications, but an interesting problem as such (though rather mathematical). To linear 
order in spin the problem is solved using an action approach, similar to a treatment of spin-^ Dirac fields 
coupled to gravity given by Kibble [6]. Further, an order-by-order construction of the canonical formalism 
with spin is given as a check. This construction is based on consistency conditions on the formalism. In 
particular it is sufficient to rely on a certain form of total linear and angular momentum expressed in terms 
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of canonical variables in order to reproduce the result of the action approach to next-to-next-to-leading 
order in the post-Newtonian approximation. The assumed form of total linear and angular momentum, i.e., 
the generators of translations and rotations, guarantees that a great part of the global Poincare algebra is 
fulfilled. The connection to the action approach is given by Noether's theorem on conserved quantities. 

Higher orders in spin correspond to quadrupole and even higher multipole corrections. Both the action 
approach and an order-by-order construction are in principle applicable to canonical formulations at higher 
orders in spin. However, only the next-to-leading order spin(l)-spin(l) level will be tackled here. This 
requires a modeling of the spin-induced quadrupole deformation, described by a single parameter for each 
object. This parameter is not only distinct for black holes and neutron stars, for the latter kind of object 
it also depends on the assumed equation of state or on other details of a particular theoretical neutron star 
model. If gravitational wave astronomy becomes available with a high enough precision in the future, one 
may hope to measure this (and maybe other) neutron star parameter. 

The results obtained here within the post-Newtonian approximation cover the next-to-leading order 
spin(l)-spin(2) and spin(l)-spin(l) conservative Hamiltonians. The conservative next-to-leading order 
spin-orbit Hamiltonian was reproduced. For maximally rotating bodies all Hamiltonians up to and in- 
cluding the third post-Newtonian order are now known. A maximal rotating body is defined to have a 
dimensionless spin (i.e., rescaled by the mass of the object and identical to the dimensionless Kerr param- 
eter for black holes) of value one, corresponding to an extremal Kerr black hole. Notice that millisecond 
pulsars (or neutron stars) and black holes can easily have dimensionless spins bigger than i (a rough ap- 
proximation for the sun yields ~ [7]). Thus spins close to maximal ones are expected to be astrophysically 
relevant. In this case the next-to-leading order spin Hamiltonians obtained here are needed for an accurate 
description of the dynamics during the inspiral phase. It was found recently in [8] that spin effects as such 
and in particular the orientations of the spins have a big impact on the event rates expected in detectors, 
especially when spins are close to maximum. 

If the fourth post-Newtonian order Hamiltonian for nonspinning objects could be obtained in the future, 
the spin Hamiltonians calculated here would be applicable to an even larger class of binaries (with smaller 
spins). Notice that the effective one-body approach for nonspinning objects, see, e.g., [9], is able to cover 
such higher post-Newtonian orders by calibration to numerical relativity and further provides predictions 
for the full waveforms, including merger and ringdown phases. An extension of the effective one-body ap- 
proach to spinning objects is possible [10]. Subsequent implementation of higher order spin Hamiltonians 
seems to be interesting, and was already performed for the next-to-leading order spin-orbit Hamiltonian 
[11]. 

Now the organization of the present thesis is given, with references to relevant published work of 
the author for certain sections (for a short review see also [12]). In section 2 spinning objects in special 
and general relativity are reviewed. Further, an overview of canonical formulations of general relativity 
is given, with emphasis on the ADM formalism and coupling to nonspinning objects. In section 3 the 
action approach to the canonical formulation of self-gravitating spinning objects to linear order in spin is 
performed [13]. An order-by-order construction based on consistency considerations is performed to next- 
to-next-to-leading order in section 4 [14, 15] as a check. In section 5 first general quadrupole corrections 
to the equations of motion and the stress-energy tensor are given [16] and then used to extend the canonical 
formalism to spin-induced quadrupole deformation at next-to-leading order [17-19]. As an application of 
the formalism, conservative Hamiltonians at next-to-leading order are derived in section 6. These are the 
spin-orbit [14] (derived earlier by Damour, Jaranowski, and Schafer), spin(l)-spin(2) [20], and spin(l)- 
spin(l) Hamiltonians, the latter was first derived for black holes [17, 18] and later for compact objects in 
general (including neutron stars) [19]. Finally, conclusions and outlook are given in section 7. 

Lower case Latin indices from the beginning of the alphabet (a, b, . ..) label the individual spinning 
objects and then consequently take on values from one to the number of objects. Three different frames 
are utilized in this thesis, denoted by different indices. Greek indices (a, fi, ...) refer to the coordinate 
frame, upper case Latin indices from the middle of the alphabet (/, J, . . . ) belong to a local Lorentz 
frame, and upper case Latin indices from the beginning of the alphabet (A, B, . . . ) denote the so called 
body-fixed Lorentz frame. Lower case Latin indices from the middle of the alphabet (i, j, . . . ) are used 
for the spatial part of the mentioned frames and are running through i = 1, 2, 3. In order to distinguish 
the three frames when splitting them into spatial and time part, we write a = (0), (i) for Lorentz indices 
(or a — (0), (1), (2), (3) in more detail), A = [0], [i] for the body-fixed frame, and /i = 0, i for the 
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coordinate frame. Indices appearing twice in a product are implicitly summed over its index range, except 
for label indices of the objects. Round and square brackets are also used for index symmetrization and 
antisymmetrization, respectively, e.g., = ^{A^ + A v ^). Partial derivatives are denoted by 9 M or 

by a comma as an index Similarly, the 4-dimensional covariant derivative is written as iu and the 
induced 3-dimensional one as A 3-dimensional vector is also written in boldface, e.g., x. The signature 
of spacetime is taken to be +2. Units are such that the speed of light c and the gravitational constant G 
are equal to one. Other symbols are defined in this thesis on their first occurrence. For convenience also a 
summary of defined symbols is given in the appendix. 

2 Preliminaries 

This section gives a short review of the achievements regarding spin in the theory of relativity as well as 
the canonical formulation of general relativity. Emphasis is put on the problems to be solved if one aims at 
a canonical formulation of self-gravitating spinning objects in the pole-dipole approximation. 

2.1 Spin in Special Relativity 

Spin already has very interesting properties in special relativity. Its canonical structure is obtained here 
as a consequence of the Poincare algebra by introducing the spin as a specific part of the total angular 
momentum. 

2.1.1 Center, Spin, and Mass Dipole 

The 4-dimensional total linear momentum P M and total angular momentum J**" = —J"^ of a physical 
system are conserved quantities due to Poincare invariance. The 4-dimensional total spin tensor S^ v can 
then be defined by 

jnu = Z n pv _ pl i z u + s »v (21) 

That is, spin is the difference of total angular momentum and its orbital part. However, a different choice 
for the yet arbitrary center Z M of the system will result in a different spin S^" (with J**" being unchanged). 
This just expresses the dependence of angular momenta on the choice of a reference point. Separating time 
and space components 

./" XT' I''/' ■ S" . J i0 ^Z^-PH + S 10 , (2.2) 

one infers that the spin transforms as 

gij S H + sz ipj _ P i 5z j ^ S i0 gio + sz i E ^ (2 3) 

under a change of the center Z l — > Z l — SZ l . E = P a is the total energy and t = Z° the time coordinate. 
Notice that J M is the total mass dipole of the system at t — relative to the coordinate origin, so (2.2) tells 
us that S l ° is the mass dipole relative to the center Z l . This explains the transformation property (2.3). 
One may also describe the 3-dimensional spin as the flow dipole and S^ u as the 4-dimensional dipole 
moment of the system relative to the center Z % . 

By its definition (2.1), transforms as a tensor under Lorentz boosts, with interesting consequences. 
In classical mechanics the center of mass, i.e., the center for which the mass dipole vanishes, is independent 
of the reference frame. In special relativity such a center can in general not be found. Under a Lorentz 
boost all components of S^ 1 " transform, so if the mass dipole S l ° vanishes in one reference frame, it will 
only be zero in all others if the system has no spin, = 0. A nice graphic interpretation is given by figure 
1. Notice that a spinning system in special relativity has a minimal extension of the order S/M orthogonal 
to the axis of rotation [21, 22]. Here S is the spin length, 2S 2 = S^ u , and M is the rest mass of the 
system, M 2 = —P^P^. In general relativity, S/M is the radius coordinate of the ring singularity of Kerr 
spacetime [23]. 
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fast & heavy 



Figure 1 : If a spinning spherical symmetric 
object moves with a velocity V 1 to the 
left, its upper hemisphere moves faster 
with respect to the reference system than 
its lower hemisphere. Thus the upper 
hemisphere has a higher relativistic mass 
than the lower one — the object acquires 
a mass dipole PAZ*. [22] 




slow & light 



However, by virtue of (2.3) one can always choose the center Z' L such that the mass dipole S M vanishes 
in one specific reference frame characterized by a timelike vector f^. That is, the center is then the center 
of mass as observed in this frame. It holds 

S^U = , (2.4) 

which is the so called spin supplementary condition. This condition fixes the center and ensures that the 
spin tensor has three independent components only. Basically three important such conditions can be 
found in the literature [22, 24], 

U=P„, or S> a 'P v = Q, (2.5) 

U = -%, or §i* = 0, (2.6) 

f M =P ll -MS%, or £""P„ - MS»° = . (2.7) 

In the following, we will indicate center and spin belonging to the second condition [21, 24] by a tilde, 
Z % and S pv ', a hat relates to the third condition [24, 25], Z % and S pv ', while center and spin of the first 
condition [26] are just denoted by Z % and S pu . We call Z % the center of inertia, Z % the center of mass, and 
Z % the center of spin [22]. Notice that the first condition is manifestly covariant, and is called co variant 
spin supplementary condition here. A different covariant condition is discussed in section 2.2. The third 
condition is called canonical spin supplementary condition, which will be explained in the following. 



2.1.2 Poincare Algebra 

The Poincare group is one of the most important groups in physics. Its generators P p and J 1 *" obey the 
Poisson bracket realization of the well-known Poincare algebra 

{P p ,P»} = 0, {P",J pa } = -Ti pp P a + r} p ' a P p , (2.8) 
{ J pu , J pa } = -rj up J pa + rj pp r a + rj afl J pv - rf v J pp , (2.9) 

where i] pLI is the Minkowski metric. Splitting space and time one gets, see, e.g, [27], 

{Pi,Pj}=0, {Pt,E} = 0, {J,,P} = 0, {G i ,P j }=E5 ij , {Gi,E} = Pi, (2.10) 
{Jj, Pj} = eijkPk, {Ji, Jj} = tijkJk, {Ji,Gj} = eijkGk, {Gi,Gj} = —eijkJk, (2.11) 

with the total angular momentum vector Jj = ^eijkJ-* and the 3-dimensional Levi-Civita symbol e^-fc. 
The boost vector J l ° has an explicit dependence on time t, which was split off as 

J m = G l -PH. (2.12) 

This defines the vector G i , which is related to the spin supplementary condition S p0 = with center of 
mass Z l by G i = Z l E, cf. (2.2). 

Notice that in general relativity total linear and angular momentum can be defined for asymptotically 
flat spacetimes as global quantities by certain surface integrals. In this case all considerations of this and 
the following section remain valid in full general relativity, see section 4.1.3. 
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2.1.3 Canonical Structure 

Using Z l = G l /E, Sy = J l0 - Z l Pj + PiZ j and the Poincare algebra (2. 10, 2.11), the Poisson brackets 
between Pi,Z % , and Sij follow as 



{Z*,P j } = 8 ij , {2T,Z>} = -j*L, {^,Z fc } = -^ + ^, (2.13) 

{Sij,Ski} = VkiSji — VkjSu — VuSjk + VijSik , (2.14) 

all other zero, where 

P P P P 

and tVy is the Kronecker symbol. 

Now we proceed to the canonical spin supplementary condition (2.7), which can be written as (E + 

M)S i0 = S ij P r From (2.3) and S i0 = we get 

Z l - Z l = 8Z l = = _*_!». (2.16) 

£ E(E + M) v ; 

Having SZ\ Eq. (2.3) relates S lj and S« by 

o _ c 1 P]PkSki PjPkSkj 

*ij ~ + E{E + M) E{E + M) ■ 
Contraction with _P; leads to EPiSij = Ml', S, , . Finally, in terms of SV, one has 



+ Af(£ + M)' J ~~ ^ M (E 1 + M) M(E + M) ' 1 



The Poisson brackets (2.13, 2.14) transform into 1 

{Z l ,Pj} = 5ij , {Sij, Ski} — SikSji — SjkSu — SuSjk + SjiSik , (2.19) 

all other zero. Thus Z % , Pj, and S!y are canonical variables. This realization is due to Pryce [24, 28]. 
Newton and Wigner further showed that Z % is the only center with this property [25]. 
Similarly, we can proceed to the covariant spin supplementary condition (2.5) by 



7, _ yi _ PkSki „ _ A , PjPkSkj Pj Pk Ski 

+ M(E + M)' ~~ j M(E + M) M(E + M) ' 1 



and find the Poisson brackets, 

{Z\Zi} = V ik V jl ^, {S iJ ,Z k } = '^^(P m S n k+PnSk m ), (2.21) 
{Z\ P 3 } = 8 l3 , {S&Su} = V£S iX - PrfSu - V^Sjk + V^Siu , (2.22) 



all other zero. 

To conclude, there are several possibilities for spin supplementary conditions and centers, however, 
only (2.7) leads to canonical variables, (2.19). This is an important fact for a canonical formulation of spin 
in general relativity. 

2.2 Spin in General Relativity 

It is well-known that spin in general relativity leads to certain gravitomagnetic effects, see, e.g., [29]. In this 
section the pole-dipole approximation for compact objects is introduced, providing an analytic description 
of spin in general relativity. 



Notice that Poisson brackets with M were calculated according to its definition M = E — Pj Pj. 
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2.2.1 Gravitational Skeleton 

In electrostatics, the multipole approximation of a charge density p, 



p(x) =(q- tfd, + y-D.n, -...) S(x) , (2.23 ) 



can be obtained from a Taylor series of its Fourier transform in the form 

p(k) = (q + iq% + ii 2 9 y feifcj + • • (2^)- 3 / 2 , (2.24) 

by the well-known transition formulas for the Dirac delta distribution <5(x) -H> (27r)~ 3 / 2 and partial coor- 
dinate derivative di -R> — iki. Here x = (x l ) are the spatial coordinates and k = (fcj) the corresponding 
ones in Fourier space. The quantities q, q l , and q %1 are the electric monopole, dipole, and quadrupole. The 
potential <p follows as 

<t> = -An A- V =(q- q'd, + 4/"^, - • • • J A > ( 2 - 25 ) 



2! 



x 



where A = d L di is the Laplacian and A -1 its inverse operator (with the usual boundary conditions). In 
most textbooks, the multipole approximation is derived directly for the potential or the field. Notice that 
the multipole approximation breaks down at high values of k, i.e., in the ultraviolet, or at small values of x 
in the potential. This is the reason for the divergent self-energy of the approximated charge density (2.23). 

Now the multipole approximation is applied to the stress-energy tensor T^ v . As it is desirable to have 
a manifestly covariant approximation scheme, we write 



-gT^ v = / dr 



1T6 W (^(4))|| Q + |(*"" Q %4))||(a 



i/3) 



(2.26) 



Here r is the proper time of a representative worldline z p (t), g the determinant of the 4-dimensional 
metric g^, <5( 4 ) = 5(x p — z p (r)), and are 4-dimensional covariant multipole moments. If one 
performs the t integration in (2.26) by eliminating the time part of <5( 4 ) and writes the covariant derivatives 
as partial derivatives and Christoffel symbols, then (2.26) indeed takes on the form of (2.23). Equation 
(2.26) in substance is Mathisson's gravitational skeleton [30], but in the form given by W. M. Tulczyjew 
[31]. Interestingly enough Mathisson unknowably used a test-function formulation of the delta distribution, 
years before this formulation was used by Laurent Schwartz for his mathematically rigorous Theorie des 
Distributions [32]. 

The divergent self-interactions already present in electrostatics become more severe if the field equa- 
tions are nonlinear. If the distributional stress-energy tensor (2.26) is used as a source for a nonlinear 
field equation, products of distributions will appear, which lack a mathematical definition. However, this 
problem can be overcome, as in quantum field theory, by a regularization and renormalization program. In 
particular, dimensional regularization [33] is most useful for theories involving gauge freedoms, like gen- 
eral relativity. Dimensional regularization has been employed successfully in post-Newtonian calculations 
[34-36] to a high order of nonlinearity. However, many treatments of multipole approximations in general 
relativity avoid these problems by considering (2.26) for test bodies only, which by definition are neglected 
as a source of the gravitational field. 

The relation between source multipoles related to used here and field multipoles [37] was consid- 
ered in [38]. Only for linear theories like electrostatics this relation is straightforward. 

2.2.2 Pole-Dipole Approximation 

The stress-energy tensor (2.26) must fulfill 

T^h^O. (2.27) 

This corresponds to Mathisson's variational equations of mechanics [30] and imposes certain conditions 
on the multipole moments. In the pole-dipole approximation only monopole t^ v and dipole t pva are kept 
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in (2.26). Evaluating (2.27) one sees that t^ v and t^ va can be expressed in terms of a vector p^ and an 
antisymmetric tensor S^, which have to fulfill the dynamic equations 



dr 



(2.28) 



with u M = 4f-, D the 4-dimensional covariant differential, and R, a the 4-dimensional Riemann tensor 



defined by 

for an arbitrary a M . The stress-energy tensor can be written as 



f—g T ^ v = \ dr 



(2.29) 



(2.30) 



p^ and 5^" are the linear momentum and spin of the object and now play the role of monopole and dipole 
moment. Their equations of motion were already derived by Mathisson [30] within his manifestly covariant 
formalism, albeit restricted to a specific spin supplementary condition. In the general form (2.28) they were 
first given by Papapetrou [39], however, his method was not manifestly covariant. W. M. Tulczyjew gave a 
derivation of (2.28) as well as of the stress-energy tensor (2.30) in a manifestly covariant way [31], using 
essentially Mathisson's method. Further important rederivations have been performed in [40, 41]. Higher 
multipole corrections will be discussed in section 5.1. 

Obviously a spinning object in general relativity does not follow a geodesic. For test bodies this effect 
can be studied numerically, see, e.g., [42, 43]. Further, without giving a relation between p p and u M , the 
system of equations (2.28) is not closed. 



2.2.3 Spin Supplementary Condition 

A spin supplementary condition (2.4) must be preserved in time. Using (2.28) this leads to a relation 
between p^ and u M [44], 

p" = — !— (-f^u" + S^^fr-) , (2.31) 



~f a u<* V dr 

and thus, for a suitable f v , closes the system of equations (2.28). A good spin supplementary condition is 
the covariant one, 



S^pu = , (2.32) 

or = pp, which has been suggested in the context of general relativity in [31]. Indeed, this condi- 
tion guarantees existence and uniqueness of a corresponding worldline z p (t) [45]. The mass quantity m, 
p M p M = — m 2 , and the spin length S, 2S 2 — S^S^, are conserved for this condition. A covariant con- 
dition has the advantage that the relation between p p and m m (2.31) is manifestly covariant. However, also 
the noncovariant condition S^ = was applied in general relativity [46]. 
A different covariant condition is given by 

S tlv u v = Q, (2.33) 



or ffj_ — u^, which was used in both special [47] and general relativity [30, 48]. While there are no serious 
objections to use this condition, as it closes the system of equations (2.28), it has some features which are 
usually not wanted. The condition (2.33) does not uniquely specify a worldline. Instead, the worldline 
depends on the choice of initial conditions and in general performs a kind of classical Zitterbewegung 
around the worldline defined by (2.32), see [31, 43]. As quadrupole corrections are needed to describe a 
black hole at the quadratic level in spin [37], we will only consider the pole-dipole approximation at linear 
order in spin here. Then the conditions (2.32) and (2.33) are fully equivalent and it holds p^ = mu^. 
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As a generalization of the canonical spin supplementary condition (2.7) to general relativity one could 
take 

S^pu + mS^riu = , (2.34) 

with some timelike unit vector n„ . However, it needs to be proven if or under which conditions (2.34) leads 
to canonical variables. 

Finally, contraction of the first relation in (2.28) with u v leads to the well-known formula 

p» = -u vV »u» - -X-^u v . (2.35) 
dr 

This relation, however, does not close the system of equations (2.28), it just is a component of (2.28). 

2.3 Canonical Formulation of General Relativity 

In this section the canonical formalism of ADM [49, 50] is introduced. Possible couplings to matter are 
reviewed, and point-masses are treated in detail. Finally alternatives to the ADM approach are discussed. 

2.3.1 The ADM Formalism 

The Einstein-Hilbert action of general relativity Wg is given by a spacetime integral over the Lagrangian 
density C G as 

H^cW = / d 4 x£ G , C G = ^V~9R (A) , (2.36) 

where is the 4-dimensional Ricci scalar. Alternatively the action can be varied with respect to the 
tetrad field e/ M instead of g^ v , see section 3.1.3. In order to find a canonical form of this action it is 
convenient to perform a splitting of spacetime into a stack of 3-dimensional hypersurfaces with constant 
time coordinate t. In these coordinates the unit normal vector n M , n^n^ = —1, of the hypersurfaces has 
the components 

1 

N 

where N is the lapse function and N l the shift vector. With the help of the projector 2 



n M = (-N, 0,0,0), or = -(1, -N*) , (2.37) 



^ =g ^ + n ^ = ( \ , ( 2.38) 

this splitting can be constructed in a geometrical way, see, e.g., [51]. The 3-dimensional hypersurfaces 
have an induced metric = jij, with jikj k ^ — Sij, a Riemann tensor Rijki, a Ricci tensor FUj, and a 
Ricci scalar R. These quantities are intrinsic geometric objects of the hypersurfaces, whereas the extrinsic 
curvature 

Ka = -n m) = (- 7 y, + 2N k ; {llj)k ) , (2.39) 

depends on their embedding in spacetime. 

Applying this splitting of spacetime to the Lagrangian density Cq leads to 

Cg = J^ N Vl [ R + K v K * J - (Hj Kij ) 2 ] + (td) , (2.40) 

where (td) denotes a total divergence, which is neglected for now. Instead of varying with respect to the 
ten independent components of g M „, we now use 7^ , N, and N l . Notice that no time derivatives of N and 
N l appear. In order to obtain a canonical formulation we have to introduce the field momentum 

tt« = 167T— = VtXtV ~ f k l 3l )K kl , (2.41) 



2 Notice that = n M 7 , " / = —N"j 0l/ and thus 7 " = for our choice of the time coordinate. 
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where (2.39) was used. This can be inverted as 
1 



Kij = 



2^7 



!,:! 



The Legendre transformed Lagrangian density then reads 



1 



Cg = -r 7r«7 iJl0 - NH + N l nf ia + (td) , 

1D7T 



field 



field 



1 

8tt 



(2.42) 

(2.43) 
(2.44) 



and the action is additionally varied with respect to tt 2 - 7 now. Notice that N and N % play the role of 
Lagrange multipliers after Legendre transformation, the corresponding constraints are the vanishing of 
"H field and Uf ld . 

A subsequent gauge fixing is subtle as it requires a fine-tuning of the action, see, e.g., [52]. As shown 
in [50, 53, 54] by different methods, see also [55], one must replace the total divergence in (2.43) by 
— for asymptotically fiat spacetimes, where £$ = — "fjj.i- This is related to the total energy E 

of asymptotically flat spacetimes by 



E 



1 

167T 



(2.45) 



where f d 2 s; denotes an integral over the asymptotic boundary of a spatial hypersurfaces at fixed time. This 
ADM energy will turn out to be the generator of time evolution after gauge fixing. For further discussion 
of boundary terms in the action of general relativity, also for the case of not asymptotically flat spacetimes, 
see, e.g., [56]. However, for asymptotically flat spacetimes the gravitational Hamiltonian may be written 

as 



H G = / d d x {NH - N l nr d ) + E[ ll3 ] . 



(2.46) 



Indeed, the action has the canonical structure momentum 7r u times velocity 7^ minus Hamiltonian Hq. 
Variation thus results in Hamilton's equations 



— — = -16-7T- = {vr 3 , H G \ . 

at djij 



9j: 



SHr 



(2.47) 



where S denotes the variational derivative here and the equal-time Poisson brackets are given by 

{7y(x),7r fc V)} = 16»r«k«fyi*(x - x') • (2.48) 

Before gauge fixing, the surface term E has no impact on these field equations, which could be obtained 
from local variations 3 . As further explained in section 2.3.3, the gauge fixing is accompanied with solving 
the constraints % field = and H^ eld = 0, so Hq then turns into the ADM energy E. To make this more 
concrete, we choose the ADM transverse-traceless gauge conditions 



d o (Hi - \lkk$ij ) = , tt" = , 
in which the transverse-traceless decomposition of 7^ and tt 1 ^ may be written as 



Hi = [ l 

nr<3 - Zti 



r ij"TT 



(2.49) 

(2.50) 
(2.51) 



3 However, one should not constrain to local variations for asymptotically flat spacetimes [54]. 
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where hJJ and 7r uTT are transverse-traceless, e.g, hJJ = hjj, = 0, and the longitudinal tt 1 ^ is related to a 
vector potential fr* = A -1 ^ 1 ^ by 

tt« = r j + fr^ - ^%^ fc , fc - -A" 1 ^,^ . (2.52) 

The advantage of this gauge is that in (2.50) there is a trace term but no longitudinal part related to a vector 
potential, while in (2.51) it is the other way around. Because of the orthogonality of the individual parts of 
the transverse-traceless decomposition, the kinetic term 71^7^0 in the action turns into 7r uTT /iJJ Q . Then 
only the transverse-traceless parts remain dynamical variables. Now the four field constraints can be solved 
for the four nondynamical variables cf) and tt 1 in terms of hJJ and 7r I:,TT . An analytic solution for <\> and fr*, 
however, can in general only be given in some approximation scheme. Notice that ADM introduced two 
slightly different gauges [50], the one used here was actually seldom used by ADM themselfes. However, 
the gauge used here is better for applications, as the form of the trace term in (2.50) is adapted to the 
Schwarzschild metric in isotropic coordinates (with obvious advantages for perturbative expansions). The 
action turns into 

W G [hJ?,% ijTT ] = / d^xTr^hJJv - { dtH ADM , (2.53) 



where the ADM Hamiltonian Hadm is just the ADM energy E expressed in terms of the gauge-reduced 
canonical variables hJJ and ir l: > TT , 

H ADM = E[h]J, ^ TT ] = -^-Jd 3 x Acf>[hJJ, ^ TT ] . (2.54) 

Notice that the surface integral (2.45) was written as a volume integral now and the asymptotic behavior of 
(j> was used. The action must be varied only with respect to the independent components of hJJ and tt^ tt , 
which is ensured with the help of the transverse-traceless projector 



(2.55) 



-(Su-A^dMdij-A^didj)}. 
The Poisson brackets after gauge fixing correspondingly read 

{^/(x), tt wtt (x')} = 16^5™<5(x _ x / } . (2 56) 

2.3.2 Matter Couplings 

Point-masses are the simplest kind of matter that can be coupled to general relativity. Its contribution to 
the action is just 

W M [9„u,z^] = J drL M , L M = -myj ' -g^z^u" . (2.57) 

This action is invariant under a change of the parameter t, which simplifies the variation as no constraint 
of the form u^u^ — — 1 is needed, m is assumed to be a constant. Variation of the action leads to the 
equations of motion 



D 

d7 



ur 



= 0. (2.58) 



These equations only have a unique solution if a gauge for r is chosen. The Einstein field equations now 
have a source ', 



1 „„„ _ „..„ ... „,..„ JW, 



Rft - §<?^%) = SttT^ , with J—gT^ = 2^ , (2.59) 
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and R% the 4-dimensional Ricci tensor. The singular stress-energy tensor density reads explicitly 



V=9T^ = / dr S {4) . (2.60) 

The 4-dimensional momentum is introduced as 

It obviously holds 4 

a r 

Im = -^V = ■ (2.62) 

Thus a Legendre transformation leads to a vanishing canonical (i.e., defined as usual) Hamiltonian. Its 
place is taken by the mass-shell constraint 

p^ + m 2 = , (2.63) 

which has to be added to the action via a Lagrange multiplier A(r), as further explained in the next section. 
This constraint is a consequence of the inability to express uniquely in terms of p M , which in turn is due 
to invariance under reparametrization, or gauging, of r. Indeed, it is a common feature of reparametrization 
invariant actions that the canonical Hamiltonian vanishes and the time evolution is instead generated by 
certain constraints. As seen in the last section this also holds for general relativity, whose action is invariant 
under reparametrizations of spacetime, or general coordinate transformations. 
Up to now the matter action was transformed into 5 

W M [9»», z^ l ,p„, A] = J dr (p^u" - H M r) , H Mt = Hg^P^Pu + m 2 ) . (2.64) 

Notice that the Hamiltonian Hm t generates an evolution with respect to the arbitrary parameter r. Fur- 
ther the variation Sp p leads to u M = 2Xp^ and from the mass-shell constraint (2.63) one thus has A = 
■^yf—u^u^. It may be checked that the equations of motion forp M and the stress-energy tensor are equiv- 
alent to the ones above, which justifies the Legendre transformation in the presence of constraints. More on 
constrained Hamiltonian dynamics is discussed in the next section. By solving the constraint and applying 
the gauge choice r = z° = t, or u° = 1, the action is expressed in terms of the independent variables pi 
and z l . It holds 

Pq = (lo 11 - n-on^p^ = g iY 3 P 3 + Nnp , (2.65) 
where np = n^p^. From the constraint we get 



-n tl n l ')p^p v + m 2 = Q np = -yj m? + Y j PtPj ■ ( 2 -66) 
Further we have 

= m =n%« = ^{goi ~ N j gji ) goi=^ijN j . (2.67) 
Putting all together we arrive at 



W M [jij ,N,N\ z\ Pl ] = Jdt { Vl z l - H M ) , H M = -Po = -Nnp - N l Pl , 



(2.68) 



4 Due to Euler's theorem, this actually holds for any Lagrangian which is a homogeneous function of degree one in the velocity 
it M . This in turn is required by reparametrization invariance. 

5 Fields within the matter action are always taken at the position from now on. 
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where a dot ' denotes the total time derivative 4?. The original Hamiltonian Hm-t vanishes by virtue of 
the constraint. Variation of the matter variables z % and pi results in Hamilton's equations with Hm as the 
matter part of the Hamiltonian. Thus z % and pi have the Poisson brackets {z l ,pj} — Sij, all other zero. As 
in the last section, the variables 7^, ir^ , N, and N l are now used for the gravitational field. 

The gauge fixing procedure is analogous to the last section, there are just certain matter corrections to 
the field constraints following from the N- and iV 4 -variations, 



= 0. 



n. 



— field _|_ -^matter q 



(2.69) 



where 



-np5 = m 2 + ^PiPj 5 , 



ft; 



(2.70) 



with 5 = <5(x — z). The first relation in (2.69) is called the Hamilton constraint, while the second one is 
the momentum constraint. The ADM Hamiltonian Hadm still results from the ADM energy by solving the 
field constraints using the gauge conditions (2.50, 2.51), but now also depends on the matter variables z l 
and pi, which have entered via source terms of the constraints. All field and matter interaction terms in the 
action, 



w - lib / d4 ^ TT ^o 



dt 



PiZ* 



(2.71) 



are contained in the ADM Hamiltonian, or the ADM energy. This is a unique feature of general relativity, 
and still holds for couplings to other matter and even other fields [53]. 

Finally we review the most important couplings of matter and fields to gravity that have received a 
canonical formulation, see also [57]. Besides for point-masses [50, 58], such canonical formulations were 
found for fluids [59], massive scalar fields [55, 60], spin-^ Dirac fields [6, 61-64], and gauge spin-1 
fields, including Maxwell [60, 65] and Yang-Mills [66]. Problematic from a canonical point of view are 
derivative-coupled theories [57], like Dirac fields and also pole-dipole objects. It is thus fortunate that the 
sought-for canonical formulation of pole-dipole objects will be seen to resemble to Dirac fields coupled 
to gravity, for which a canonical formulation was found. Though the classical spin of pole-dipole objects 
is not restricted in its size, we consider pole-dipole objects only at linear order in spin here. This means 
the spin is treated as an infinitesimal quantity and thus formally takes on the smallest (nonzero) classical 
value, which seems to give rise to similarities to the minimal (nonzero) quantized spin -| of Dirac fields. 
Thus the achievements on canonical formulations of Dirac fields coupled to gravity served as a very useful 
guide here, in particular the paper of Kibble [6]. However, an additional problem to be solved for spinning 
objects in general relativity concerns the canonical spin supplementary condition. So far the canonical 
formulation of spinning objects was found for test-bodies in an external gravitational field [67], see also 
the very recent work in [44] . 



2.3.3 Other Formalisms and Constrained Hamiltonian Dynamics 

Before gauge fixing, general relativity possesses a canonical formulation in the presence of the constraints 
H = and Hi = 0. There exists a general framework to handle such a constrained Hamiltonian dynamics, 
which was developed most notably by Dirac as a general route to canonical quantization [68], see also 
[52, 69, 70]. Further important work was done by Bergmann and his collaborators, but focused on general 
relativity and its canonical quantization [71]. Though Dirac also considered the canonical formulation of 
general relativity [72], his approach is formulated in a very general way. Early work on this subject can 
even be traced back to Rosenfeld [73]; for a historical review, see, e.g, [74, 75]. A particular important 
achievement of ADM for canonical general relativity was the identification of the ADM Energy as the 
Hamiltonian after gauge fixing [49, 76]. Yet another canonical treatment of general relativity was given 
by Schwinger [55]. This formulation is similar to the ADM one, essentially only different variables were 
used and many more such reformulations are possible. A further very appealing formulation was given by 
Ashtekar [77], in whose variables the gravitational constraints considerably simplify, and which forms the 
basis of loop quantum gravity, see, e.g., [78]. 



2.3 Canonical Formulation of General Relativity 



15 



We will now summarize some of the results of Rosenfeld, Dirac, and Bergmann on constrained Hamil- 
tonian dynamics. In the last section the mass-shell constraint (2.63) manifests the inability to uniquely 
express the velocity in terms of the corresponding momentum p^. The standard route to a Hamiltonian 
seems to be impassable in such a situation as the Legendre transformation can not be applied in its usual 
way. The solution, however, is simple. The Legendre transformation may formally be performed as usual 
if one adds the emerging constraints via Lagrange multipliers to the action. The additional degrees of free- 
dom introduced by these multipliers correctly parametrize the ambiguity present in the relation between 
velocities and momenta. Further, it can be shown that the dynamics of the transformed action is equivalent 
to the dynamics of the original action. The constraints arising at this stage are entitled as primary and the 
Hamiltonian is called the total Hamiltonian (or Dirac Hamiltonian), as it includes the primary constraints 
via Legendre multipliers. 

The next step in the analysis of constrained Hamiltonian dynamics consists of evaluating the consis- 
tency requirement that all primary constraints must be preserved under the time evolution given by the 
total Hamiltonian. Of course, some of the resulting consistency conditions can be identically fulfilled or 
lead to contradictions (then the dynamics must be considered as inconsistent). Moreover some conditions 
are restrictions for the Lagrange multipliers appearing in the total Hamiltonian. Due to linearity of the 
total Hamiltonian in the Lagrange multipliers, these restrictions are actually linear equations. Further, one 
might also obtain new (independent) constraints from the consistency conditions. Such new constraints are 
called secondary constraints. For these new constraints the same consistency requirement applies, and one 
is eventually lead to further conditions on the Lagrange multipliers and/or to further secondary constraints 
and so on. Finally, one ends up with a complete set of constraints and linear equations for the Lagrange 
multipliers. 

The linear equations for the Lagrange multipliers can be used to eliminate certain linear combinations 
of these multipliers from the equations of motion. The usual situation known from courses on classical 
mechanics is that all multipliers are uniquely fixed. However, in the general case some combinations of 
Lagrange multipliers could remain unfixed and thus remain as arbitrary degrees of freedom in the equations 
of motion. The interpretation is that these degrees of freedom are physically irrelevant and correspond to a 
gauge freedom of the theory. That is, the corresponding independent Lagrange multipliers can be chosen at 
will, interpreted as choosing a gauge. Hamiltonian formulations of gauge theories will inevitably involve 
constraints. 

The Lagrange multipliers enter the total Hamiltonian together with the primary constraints. Instead 
of characterizing the gauge freedom of a theory by undetermined combinations of Lagrange multipliers, 
one can give a description in terms of corresponding primary constraints. For this purpose it is useful to 
introduce the notion of first class and second class constraints. First class constraints are defined to have 
vanishing Poisson brackets with all other constraints. A constraint that is not first class is called second 
class. In addition to being first or second class, the constraints can still be primary or secondary, and one 
thus has four categories of constraints now. An important fact is that the number of independent primary 
first class constraints is equal to the number of unfixed Lagrange multipliers in the equations of motion and 
thus to the number of gauge degrees of freedom. 

Not only the primary first class constraints but also all secondary first class constraints are related to 
gauge symmetries [79, 80] (at least under certain reasonable conditions), see also [75]. To be more precise, 
all first class constraints, primary as well as secondary, appear in the generators of gauge symmetries on 
phase space. The algebra of first class constraints is therefore related to the algebra of gauge symmetry 
generators of the theory. For general relativity, the algebra of first class constraints reads [53, 55] (at least 
for the vacuum case and for coupling to point-masses) 



If one goes to the constraint surface by % = = %i, then the right-hand sides vanish. Thus % and %i are 
indeed first class. In order to relate this algebra to 4-dimensional diffeomorphism invariance one should 
include lapse and shift as well as corresponding momenta into phase space [79]. It should be noted that 
though the total Hamiltonian of general relativity (2.46) is composed of the first class constraints and looks 



{%(x),H(x')} 
{^(x),H,(x')} 



[%(x) 7 «(x) + H 4 (x') 7 4 V)] djS(x - x') , 
W(x)ftJ(x-x') ) 

Hj(x) <9 2 <5(x - x') - Wi(x') dj5(x - x') . 



(2.72) 
(2.73) 
(2.74) 
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quite similar to the generator of gauge transformations, the time evolution given by this Hamiltonian is not 
just a gauge effect, see, e.g, [81]. 

One can elaborate more on the distinction between first and second class constraints. Obviously one 
can recombine the whole set of constraints into some equivalent set. We consider the case that such a 
recombination brings as many constraints as possible from the second class into the first class. One can 
then show by a reductio ad absurdum that the matrix c a \, — {ip a ,ipb}, where ifj a are the constraints that 
remain second class after recombination 6 , is invertible, det(c a b) 7^ 0. The Dirac bracket {A, B}* between 
two phase space functions A and B is then defined by 



This bracket satisfies the laws known from the Poisson bracket. Further, it leads to the correct equations 
of motion together with the total Hamiltonian. The Dirac bracket can thus be used as a substitute for 
the Poisson bracket. However, whereas one may use the constraints only after all Poisson brackets were 
calculated 7 , the second class constraints ip a = can be used before an application of the Dirac bracket 
without changing the result (e.g., one has {A, ip a }* — for all A and ip a ). If one restricts to use the Dirac 
bracket instead of the Poisson bracket, one can use the second class constraints tp a = to solve for certain 
phase space variables and eliminate them from all quantities. Then one has performed an actual reduction 
of the degrees of freedom. Within this formalism, gauge conditions are constraints added by hand that 
bring all (or just some) first class constraints into the second class. The reduction of degrees of freedom 
via gauge fixing then follows with the help of the Dirac bracket in a straightforward way. 

3 Action Approach 

In this section an extension of the ADM formalism for point-masses to the pole-dipole approximation is 
obtained linear in spin. The derivation is based on a corresponding extension of the point-mass action. 

3.1 Action of the Spherical Top 

It is remarkable that equations of motion (2.28) and stress-energy tensor (2.30) in the pole-dipole approxi- 
mation are independent of the specific object, i.e., are the same for black holes and neutron stars. It is thus 
expected that any specific action for a spinning object coupled to general relativity will contain (2.28) and 
(2.30) to some approximation (e.g., linear in spin). The pole-dipole action found here will be based on the 
simplest spinning object imaginable — the spherical top. 

3.1.1 Newtonian Case 

The spherical top is well-known in classical mechanics. However, we review it here in a way that allows 
an easy transition to the special relativistic treatment in [69, 82]. We consider in this section a top with its 
center of mass resting at the coordinate origin. The center of mass motion can be added easily. The top 
can be described as a rigid body consisting of many point-masses labeled by an index a, with positions 
z % a and masses m a . In terms of body-fixed (constant) coordinates Za it holds z l a (t) — A^ l (t)za\ with 
a time-dependent rotation matrix A^ k ' l A™ — Sij. Here and in the following we will indicate 3- 
dimensional indices in the body-fixed coordinate system by square brackets. The rotation matrix can be 
expressed in terms of three independent angle variables, A^ = A^(ipi, ip2, ^3), e-g-, the Euler angles. 
The antisymmetric 8 angular velocity tensor is given by Q % i = A^'A™. A spherical top is completely 
characterized by one moment of inertia /, it holds 2 ^ a m a z$Za — I Sij. The Lagrangian of the free 
spherical top then reads 



{A, B}* = {A, B} - {A, fayic-^abtyb, B} . 



(2.75) 




(3.1) 



a 



The indices a and b label constraints in this section. 
7 For a more detailed exposition it is useful to introduce the concepts of weak and strong equality. 
8 The antisymmetry immediately follows from the time derivative of A.W*A.Wi = 8ij . 
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where ST' = ieyfcfP is the usual angular velocity vector. The spin is the generalized momentum of the 
angular velocities of the form 



dL 

SiJ = = ISly ■ (3.2) 



Legendre transformation leads to 



L = -S ij Q ij -H, (3.3) 



with the Hamiltonian H(A^ , Sij) — jj SijSij. This specific Hamiltonian is actually independent of AW- 7 . 

In order to derive the general Euler-Lagrange equations of the Lagrangian (3.3) we are not varying the 
independent angle variables, but instead use 58 l i = A^SA™ as independent variations. Notice that 59^ 
is antisymmetric, and thus indeed corresponds to three independent variations of the angle variables. The 
result is 

- AW * A[feb ' - 2 |f- • ^ = 2 %%I* ~ AW'^j + AWi^j . (3.4) 

These are Hamilton's equations for and Sij. The Poisson brackets fulfill 
8 A 

A = {A,H} + — , (3.5) 
for a general quantity A. Comparing with (3.4) we can read off 

{AWaM<} = 0, {AW,S k i} = A® k 5 lj -AW5 kj , (3.6) 
{Sij, Ski} — SikSji — SjkSii — SuSjk + SjiSik ■ (3.7) 

Alternatively one could use canonical variables based on the angle variables, 



[r,pi\ = ^r Withpf ^y*' , .3.8) 



as in most textbooks. Further, if the Hamiltonian is independent of A™', which will always be the case in 
the following, the spin length is a constant and it is possible to describe each spin by only two independent 
canonical variables instead of six contained in Sij and A^, see, e.g., [83, 84]. However, we prefer the 
variables Sij and A^- 7 here. 

3.1.2 Special Relativistic Case 

In the relativistic case there are no rigid bodies. However, one can define a top in a purely mathematical 
way [69, 82, 85] as a worldline with a Lorentz matrix Aa/j,, r/ AB Aa^Abu = Vnv, such that A^ is a pure 
rotation, 

Aau = ( ~} ° V (3.9) 



in some frame defined by f^. (Upper case Latin indices from the beginning of the alphabet refer to the 
body-fixed frame and have the values A = [0], [i].) This can be formulated as 

A mu = -., or r^ A = —^k=. (3.10) 



It holds 
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see, e.g, [69] or the next section. The spin supplementary condition belonging to (3.10) reads 

S^r = 0. (3.12) 

It will be seen in section 5.2.2 in which sense this belongs to (3.10). Notice that only three relations of 
(3.10) are independent, e.g., one could equivalently require A^/^ = only. The same holds for (3.12). 

There are many ways to implement the conditions (3.10) and (3.12) in an action approach, see, e.g., 
[47, 69, 82]. We require here that (3.10) and (3.12) are preserved under the time evolution given by the 
action and try to directly construct such an action. An alternative, rather indirect, approach would be to add 
the supplementary conditions to some action with the help of Lagrange multipliers. As well known from 
classical mechanics, this modifies the dynamics by constraint forces, which ensure that the supplementary 
conditions are preserved in time. However, one should carefully check the consistency, in particular one 
should be able to find a solution for the Lagrange multipliers. Also no further (secondary) constraints 
should appear, which would be physically unacceptable (we want to have exactly three independent rota- 
tional degrees of freedom). Finally, the Lagrange multipliers can be eliminated from the action, leading to 
a dynamics which preserves the constraints and thus to the action we try to find directly here. 



3.1.3 Minimal Coupling to Gravity 

The next logical step is a minimal coupling of the special relativistic spherical top defined in the last section 
to gravity. Such a coupling was already treated in [86] based on the developments in [82]. In [87] even 
nonminimal couplings leading to higher multipole corrections were considered 9 . Notice that [88] is not 
a further development of [87], but is a completely different action approach. More recently yet another 
approach was given in [89] with focus on an application to the post-Newtonian approximation. 
The matter variables A^ have the problem that they fulfill 



A A ^ 



(3.13) 



That is, A^ is not independent under variation of the metric. For the Dirac field, one has a similar problem 
with the gamma matrices 7^, as it holds 7^7^ + 7^7^ = 2g^ v . This problem can be overcome by writing 
A A fl = A AI ej^ and treating A AI and the tetrad field ej^ as independent variables. From 



AaiA a ., 



mj : 



(3.14) 



it is obviously now consistent that Aai and 7/ are constant under variations of the tetrad field e/ M . We have 
three bases involved, a body-fixed basis, a local Lorentz basis (denoted by upper case Latin indices from 
the middle of the alphabet), and a coordinate basis. The field equations are obtained by an unconstrained 
variation of e/ M . The metric g^ v = ej^e 1 v as well as the connection are not varied independently. For the 
variation of Aai one has to take into account the condition (3.14). 

In [86, 87, 89] matter and field degrees of freedom are not clearly separated in the action. For example, 
in [87] the equations of motion for the matter variables were obtained by adding (3.13) as a constraint 
to the action with the help of Lagrange multipliers, whereas the field equations were obtained from an 
unconstrained variation of A^. However, we need to separate matter and field degrees of freedom here, 
which is essential for the canonical reduction in the next sections. 

The covariant angular velocity in the local Lorentz basis can be defined as 



u 



A/ 



DA 



A J 



dr 



A, 



dA 



A.J 



At 



(3.15) 



Here uj i1 ij are the Ricci rotation coefficients, eiaejpu^ 1,1 = — T^^ + e K aifl eK/3, and Taj^u = \{g 



gav,fi — g^u.a) is the 4-dimensional Christoffel symbol of first kind. Notice that the covariant derivative 
does not act on indices referring to the body-fixed frame. The matter action shall be of the general form 



W M [ei^z^A AI ] = / drL M (u^n^,g^ P ))- 



(3.16) 



'This obviously goes beyond a spherical top, however, the formalism stays the same. 
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The Lagrangian Lm is restricted to depend on the velocities up and Q,^ v = ei^ej v Q IJ only, and not on the 
"coordinates" z M and A AI directly. This ensures the covariance of the action. The action shall be invariant 
under reparametrizations, so is not constrained. If we let the Lagrangian depend on the curvature tensor, 
we would include quadrupole corrections, see section 5.2. An important relation is given by (2c) in [82] or 
(9) in [87], which reads here 



du a 



dL 



A I 



2 „ 9au j 



(3.17) 



and is a consequence of Lm being a scalar 10 . Similar to section 3.1.1, the Euler-Lagrange equations are 



obtained with the help of the antisymmetric variations . 



A/<5A 



A3 



e.g.. 



DS9 



u 



i j 



(3.18) 



The S9 IJ -variation then leads to 



D 
d7 



dL 



A I 



dL 



A I 



(3.19) 



The Jz^-variation is subtle as it is not manifestly covariant, see, e.g., the second last term in (3.18). This 
is due to the fact that A AI is held constant for the variation of the worldline SzP, which is not a covariant 
process (e.g., in contrast to a parallel transport of A AI to the new worldline). However, using the equations 
of motion for A AI , (3.19), and the covariance of Lm, (3.17), the result of the <5z M -variation reads 



D 



dL 



A I 



duP 



R (4)a0 v 9L M 

p v dn a ? ' 



(3.20) 



and is manifestly covariant now. Further, by virtue of (3.17) we can write (3.19) as 
' dL M ' 



dLM dLM 
-u v — — u,. 



duP 



du" 



At last, the field equations follow from the <5e/ M -variation as 
RZ - = SttT^ , with ^—gT^ = er 



A I 



5ei v 



(3.21) 



(3.22) 



where the left-hand side results from the Einstein-Hilbert part (2.36) and the stress-energy tensor density 
yJ—gT 1 * 1 ' reads explicitly 



-gT^ v = / dr 



u^g 



dL M ap 0(u v)g 



Up. 



Here the important relation (3.17) was used again and the antisymmetric part 



dr 



D f dL 



A I 



dr \dtt a f 



(4) 



vanishes, see (3.19). Indeed, (3.19) is equivalent to = 0. 

Comparing (2.28) and (2.30) with (3.20), (3.21), and (3.23) we get 



~ _ 9 dL M 



dL 



A I 



duP 



(3.23) 



(3.24) 



(3.25) 



as in the special relativistic case. It should be noted that the given derivation basically follows along the 
lines of Bailey and Israel [87], but the used variables are similar to Porto [89], which resembles to [69]. 



'"Loosely speaking, one can read (3.17) as "the number of upper indices minus the number of lower indices in Lm is zero." This 
is derived in [87] from an infinitesimal coordinate transformation. 
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However, the variables that are varied here differ from both [87] and [89]. The relation between p^ and u M 
is fixed by (3.25), which means that the action already implements a specific spin supplementary condition, 
cf. Eq. (2.31). If this would not be the case, then (3.25) should be of the form (2.35), which is impossible 
due to the assumed absence of accelerations in the action. The approach in [82, 86] includes accelerations 
of the worldline coordinate. Further, a noncovariant supplementary condition (e.g., the canonical one) will 
result in a not manifestly co variant relation between p^ and u M , (3.25), and thus one needs a not manifestly 
covariant action. An approach via Lagrange multipliers as discussed in section 3.1.2 seems to be better 
when using such conditions from the start. Here we will start with the covariant supplementary conditions 
and go over to the canonical ones later by a change of variables. 

Now we have to find a suitable reparametrization-invariant Lagrangian. An intuitive guess is (see 
section 5.2 for more elaborated considerations) 



Lm 



(3.26) 



y/-U p uP 

where and I shall be constants. Then it holds 

S M , = -^k=, p^L + kx)^, (3.27) 
and the dynamical mass m = ■J—p^p' 1 is given by 

m = m + jjS a pS a P , (3.28) 



or m = mo to linear order in spin. Then (3.27) agrees with (2.31) for / M = p^ at linear order in spin, 
which implies that the corresponding spin supplementary condition (3.12) is preserved in time. (3.10) 
only needs to be preserved to zeroth order in spin, which is also the case (see also section 5.2.2). Due 
to reparametrization invariance, Lm must be a homogeneous function of degree one in the velocities and 
Euler's theorem leads to 



LM = ^ + 1^" = P ^ + ■ (3 ' 29) 

A Legendre transformation in and W thus leads to a vanishing result. Further, the mass-shell constraint 
(2.63) follows from (3.27), but no constraint on S M „ arises from (3.27) as opposed to [69]. (Indeed, in [69] 
the action was constructed such that the constraint (3.12) arises directly from the action in this way.) Similar 
to section 2.3.2 we finally have 



WM[ei^z", P p, S^,A AJ ,A] 



dr 



1 



(3.30) 



with the function 11 Hm t containing the mass-shell constraint only, Hm t = Md^P^Pf + m2 )- This is 
the extension of (2.64) to the pole-dipole approximation at linear order in spin. We could also add the 
supplementary conditions, 



(3.31) 



to the action with the help of Lagrange multipliers. However, this will not change the dynamics as these 
(independent) conditions are already preserved in time and their Lagrange multipliers therefore vanish. 
Reference [13] immediately started with the action in the form of Eq. (3.30) without the detailed derivation 
given in this section. 



Notice that % T is not a Hamiltonian as | S^f^" in (3.30) also contains interaction terms. 



3.2 Reduction of the Matter Variables 
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3.2 Reduction of the Matter Variables 

Next a fully reduced canonical formalism is derived. For this the action is put on the constraint surface. 
That is, all supplementary conditions, constraints, and gauge conditions are solved in terms of certain truly 
independent variables that parametrize the constraint surface. The equations of motion for this reduced 
number of variables could then be obtained by varying the action with respect to these variables. However, 
we will transform the action to a new set of reduced variables such that the equations of motion can easily 
be seen to resemble to Hamilton's equations. This allows for an easy identification of the Hamiltonian and 
corresponding Poisson brackets. Thus a fully reduced canonical formalism for spinning objects coupled 
to general relativity is found [13]. Remember that the necessity for a variable transformation to obtain 
standard canonical Poisson brackets is already present in the fiat space case, see (2.20). A treatment using 
Dirac brackets (2.75) seems to be more complicated, as one has to consider the brackets for each pair of 
variables then, whereas here we only have to handle the action (a single scalar). 

The derivation sketched above is very similar to the treatment of Dirac fields coupled to gravity by 
Kibble [6]. In this section we concentrate on the matter part of the action only. 



3.2.1 Reduced Matter Action 



Similar to section 2.3.2 we solve the matter constraints (now including the supplementary conditions (3.31)) 
as 



np = n^p^ 



TiSi — XI Sui 



Pkl k] S 3l 



np 



lijnS 3 



A [j](o) = A H(.)M A [o]/ = _P_ 



,(0) ' 



III 



(3.32) 
(3.33) 



in terms of the independent variables pi, Sy, and AW'". On the constraint surface it holds Hmt = 0. 
For simplicity, we will immediately constrain ourselves to the Schwinger time gauge [55], 



e (0)n = , 



(3.34) 



see also [6, 61, 64], as lapse and shift then turn into Lagrange multipliers in the matter action [13], like in 
the ADM formalism for nonspinning objects. This gauge condition effectively reduces the tetrad e /M to a 
triad e"", it holds 



j( ) i = = ew °, 



N l = -Ne 



(o) 



M Q — N — l/e (0) ° 



;« = N j e^ 



l lJ = e( m )-e 



j > 



(3.35) 
(3.36) 



A further convenient gauge choice is r = z° = t for the yet arbitrary parameter r. In terms of the 
independent variables the matter Lagrangian (3.29) reads explicitly 



Lm — 



p l + K l3 nSi + A kl e {j)k e<j\ ti 

nS (i)PU) ~ nS U)P(i) 



2 np 



nS l . 
2np^ 



'WO') 



n p 



A [fc] WAM«) 



(3.37) 



+ A %3 e [k)i e^> jfi - J d 6 x (AH matter - N l H? ltter ) , 
with the 3-dimensional Christoffel symbols Tkji, the abbreviation A ,J defined by 



1 



tS, 



iPj 



2~' J ' 2np ' 
and the matter parts of the gravitational constraints given by 



= -npS- K^^^S - {nS k S). k , 
np 



(3.38) 



(3.39) 
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(ft + K l3 nS>)5 + ^ k S lk S + l jkE ^s) 



nr ner = (ft + i^n^)<5 + \~y K S ik 5 + t""-*—=L6 ) . . 3.40! 

These coincide with the densitized projections 

^matter = ^ T ^ n » n » , ^J-tter = ^^ Tlv n v , (3.41) 

of the stress-energy tensor (2.30) at linear order in spin. For consistency this must of course be the case, 
as the gravitational constraints can also be obtained by such projections of the Einstein equations directly, 
instead of by varying the action with respect to N and N z . However, in spite of the simplifying premature 
(but only partial) gauge fixing (3.34) of the tetrad, the matter Lagrangian (3.37) is still complicated com- 
pared to the nonspinning case (2.68). In particular, the canonical structure is not immediately visible in the 
used variables. 



3.2.2 Canonical Matter Variables 

One already knows from special relativity that the variables in the covariant spin supplementary condition 
have quite complicated Poisson brackets. Thus the complicated structure of the matter action in these 
variables found in the last section is not surprising. We will now try to simplify the structure of the 
matter Lagrangian by introducing new variables, which will turn out to possess standard canonical Poisson 
brackets. These new variables are indicated by a hat. An intuitive guess from the special relativistic case 
(2.20) is 

z i = gi _ , nSi = _PJJ^ s = §i . _ + If^i. f (3 .42) 

m — np m m — np m — np 

belonging to the condition (2.34), as well as 

o(j') 



m(m — np) / 

see (3.60c) in [69]. These redefinitions replace A kl in (3.37) by the quantity A y given by 



(3.43) 



\kl _ 1 c m P(i nS 3) 



likljiA = -S i:i + p J ' . (3.44) 
2 np[m — np) 



Then the first line of (3.37) suggests to introduce a new linear momentum for the matter as 

Pi = pi + K ijn S j + i fci e (j)fc e« M - Qs fc , + MLil^I*^ ' 



(3.45) 



which reduces to pi = P i in the special relativistic case. The matter Lagrangian now turns into (still 
approximating linear in spin) 

L M = fti* + ^3 m) M m - H M , (3.46) 

where = A [k] ^A^ and 

H M = -A ij e ik)i eWj,a + J d 3 x {NH mme1 ' - N l U™ ma ) . (3.47) 

Notice that AM^') is a 3-dimensional rotation matrix, Arw^A^w) = Therefore fiWwJ is antisym- 
metric and should be interpreted as an angular velocity tensor. The action thus has the canonical structure 
momenta times velocities minus Hamiltonian Hm. The Poisson brackets for the matter part read 

{z\ Pj } = 5 tj , {A m \S m) } = AH( fe )% - A m) 5 kj , (3.48) 



3.3 Full Gauge Reduction 
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{S(i)(j)> ^(fc)C)} = $ikS(j)(i) ~~ $jkS(i)(i) — 5aS(j)(k) + fijiS(%)(k) > (3.49) 

all other zero, similar to (3.6, 3.7). It is important that all extrinsic curvature terms are eliminated from 
(3.37, 3.39, 3.40) by the redefinition of the linear momentum (3.45). Terms of this type are the reason for 
potential problems with derivative-coupled theories [57], so it is good that they disappear. This is similar 
to the Dirac field case, which can be made a nonderivative-coupled theory by a redefinition of the Dirac 
field. Further the p;-term in (3.37) was removed by the redefinition of the position (3.42). 

If we consider test spinning bodies in an external field, then one immediately gets the fully reduced 
Hamiltonian in the time gauge by inserting the metric (i.e., 7^, N, and N l ) as well as a suitable triad e^)i 
(subject only to erfAi e j — 7^ ) into (3.47). Canonical formulations of test spinning bodies were already 
obtained in [67] by a direct construction of the symplectic structure and also very recently in [44] using a 
Dirac bracket approach. In the latter paper the Hamiltonian was explicitly obtained for the Kerr metric. In 
the next section we will also be able to put the field part into canonical form. 

Given the fact that, at least in the time gauge (3.34), the supplementary condition (2.34) leads to a 
canonical spin and position variable, it seems to be simpler to immediately start with an action implement- 
ing (2.34), thus skipping the need for variable redefinitions. However, one can not be sure in advance that 
(2.34) leads to canonical variables. Further, it should be noted that only the structure of the action was 
simplified by above redefinitions. The redefinitions still have to be applied to (3.39, 3.40), making these 
expressions more complicated, see (6.33-6.35). Thus one has a conservation of trouble here and starting 
directly with (2.34) does not seem to simplify the calculation. In fact, it could be subtle to correctly imple- 
ment the noncovariant condition (2.34) into the action. However, this succeeded for test spinning objects 
in [44]. 



3.3 Full Gauge Reduction 



The discussion of the field part is not as simple as for nonspinning objects. First, we need the tetrad form 
of the ADM formalism as derived in [62]. Second, the matter action depends on the partial time derivative 
of the tetrad, which necessitates matter corrections to the canonical field momentum. Indeed, the canonical 
momentum conjugate to erjAj is given by 



8tt 



d(C G + C M ) 



de 



(*)3>0 



(3.50) 



where Cm is the density version of (3.46), obtained by introducing 5 = S(x l — z 1 ) in certain terms, and n* 3 
is still given by (2.41). Remember that (3.39, 3.40) do not contain the extrinsic curvature after redefining 
the matter variables. Legendre transformation leads to 



W = 



H = 



1 

8^ 



dt 



1 



PiZ 



:S 



d :! . 



(NH - N l H t + AyTrM) + £[7y] , 



(3.51) 
(3.52) 



where U = H md + H mMta ' and Hi = Hf eld + Hf^ 1 with (2.44) and (3.39, 3.40). In tetrad gravity one 
has the additional constraint 7r™ = 0, or 7f™ = 8ttA^S, which was added to the Hamiltonian H via a 
Lagrange multiplier = —Xji- 



3.3.1 Spatial Symmetric Gauge 

The constraint = is eliminated by a further partial gauge fixing now. The spatial symmetric gauge 
for the triad eu\j = ejj = eji is imposed, which was suggested by Kibble for a canonical formulation 
of the Dirac field coupled to gravity [6] (however, Kibble was using the Schwinger canonical formalism 
[55]). In this gauge, the triad is the symmetric matrix square-root of the positive definite induced metric, 

eijejk = Jik, or 



0«) = \/(7ij) 



(3.53) 
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Thus the triad is fully given in terms of the metric, which is now the variable to be varied. We may therefore 



define an object BM as 



?B kl -e d€mj c d6mt (3 54) 

OJkl oj k i 



which enables us to write 



1 



e (fc) 'ie(fc)j,/i = -By 7M,/i + ^7ij !jU . (3.55) 



(3.56) 
(3.57) 



The action obviously takes on the form 

w = ~ikJ dAx%ii ^° + J dt + 5%)(i)^ (i)0) - 17 

= J d 3 x (NH - iV%) + E[-y tj ] , 
with the new canonical field momentum conjugate to 7^ given by 

n ij = n ij + 8nA^d+ WnB^A^S. (3.58) 

We have thus reduced the tetrad form of the ADM formalism to its metric form, still coupled to spinning 
objects. 

3.3.2 ADM Transverse-Traceless Gauge 

Finally, the gauge fixing for the induced metric follows along the same lines as for nonspinning objects in 
section 2.3.2. We apply the gauge conditions 

9j (7ti " hkkSij ) = , tt u = . (3.59) 

However, notice that the ADM transverse-traceless condition for the canonical field momentum ft n = 
differs from the original one, tt" — 0. Correspondingly we now have the decomposition 

rt» = + 7r^ TT , #« = + ^ - i^l fe , fc - -A" 1 ^* , (3.60) 

instead of (2.51). The decomposition for the metric (2.50) is still valid. The ADM Hamiltonian then results 
from solving the field constraints H = = %i together with the gauge conditions as 

H ADM = E[z\p i J mj) ,hJJ,ir ijTI } = J d 3 a;A0, (3.61) 

and the fully reduced Poisson brackets of the field read 

{ft™(x), 7r fem (x')} = 16t4™S(x - x') , (3.62) 
all other zero. The Poisson brackets (3.48, 3.49) of course still hold. The fully reduced action finally reads 



w =TLi d4x * if " h ™ + ldi 



p, 



(3.63) 



This is the extension of the nonspinning case in (2.71). The new spin interactions enter via the ADM 
Hamiltonian H^dm after solving the constraints, which now have spin corrections in its source terms, 
(3.39, 3.40). 



4 Symmetry Generator Approach 
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As we have seen in section 2.3 and also in the last section, after all constraints as well as supplementary 
and gauge conditions have been eliminated, the Hamiltonian is given by the ADM energy depending on 
the fully reduced canonical variables. However, while it is not problematic to calculate the ADM energy 
at least to some order in a perturbative way, it will then depend on the variables appearing in the stress- 
energy tensor (2.30) and equation of motion (2.28), for which the canonical structure is not known. If 
one could someway find the transformation between these variables and fully reduced canonical variables 
with usual Poisson brackets, then the ADM energy can be expressed in terms of canonical variables and 
turns into the ADM Hamiltonian. In this section we try to construct this variable transformation order-by- 
order in some perturbation scheme by looking at certain consistency conditions. It is expected that if one 
proceeds to higher and higher orders, then one also needs to devise more and more consistency conditions. 
However, in the post-Newtonian approximation one may reach an order high enough for all currently 
relevant applications by just relying on a specific form of total linear and angular momentum expressed in 
terms of canonical variables [14, 15]. Notice that this approach is not as powerful as the action approach 
[13] discussed in the last section, however, it succeeded earlier and is still valuable at higher orders in spin 
as well as for a check of the action approach at linear order in spin. 

4.1 Symmetries and Conserved Quantities 

Now the symmetries and corresponding conserved quantities for asymptotically flat spacetimes are re- 
viewed. These conserved quantities generate their symmetries on phase space. For total linear and angular 
momentum this leads to a very specific form when expressed in terms of canonical variables. 

4.1.1 Global Rotations and Translations 

It is intuitively clear that an asymptotically flat spacetime can be transformed into a physically equivalent 
one by a 3-dimensional rotation and/or translation of each 3-dimensional hypersurface, i.e., of the whole 
spacetime. This means that asymptotically flat spacetimes posses a global symmetry 12 under rotations and 
translations, i.e., under the 3-dimensional Euclidean group. In fact, one even has a global symmetry under 
the Poincare group, which will be discussed in section 4.1.3. How the symmetry under the Euclidean group 
is represented on the coordinates crucially depends on the chosen coordinate system even in flat space. If 
the coordinate system resembles to a Cartesian one in the asymptotics, then a good guess for the symmetry 
transformation is x % — > Ayfx- 7 + a?), with x % the coordinates of the 3-dimensional hypersurfaces, a 1 a 
constant vector describing a translation, and a rotation matrix Ay . The rotation matrix is parametrized by a 
constant antisymmetric matrix cj ij = — uj^ 1 in the form A = e w . A field, e.g., the induced metric 7^, then 
transforms as 

7y-(x) -» AikAjakiiA^x - a) , (4.1) 

where the vector a has components a 1 . However, for this transformation to be a global symmetry and not 
just a particular gauge transformation, the gauge conditions must be invariant under this representation of 
the Euclidean group. This is indeed fulfilled for the ADM gauge conditions (2.49) or (3.59) (remember 
that a 1 and Ay are constant). Further, the local basis shall rotate the same way as the coordinate basis, 
i.e., e^jfx) — > AifcAj/e(fc);(A -1 x — a). We assume here that the tetrad was reduced to a triad with the 
help of the time gauge (3.34), as in the action approach. The triad gauge shall be compatible with this 
transformation property, which is the case for (3.53). 

Looking at infinitesimal transformations, i.e., a 1 and u lJ shall be small, it holds 

x l ->■ x l +a l + w*V , (4.2) 

or for a tensor field (4.1) 

Hj -> Ho ~ a k d k y l3 - w kl x l d k 'Yij + uj tk ^ k] + w J ' fe 7 ifc . (4.3) 

12 A global symmetry depends on parameters which may not vary over spacetime. 
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This is just the Lie-shift given by the infinitesimal coordinate transformation (4.2), i.e., — > 7y — 
£&x k lij- Similarly, the canonical variables transform as 

z\ -> z l a + a 1 + u 11 z{ , p ai -> p ai + uj' ij p aj , (4.4) 
AWO) AJP + w»' fc AWW , ,§ o( . )W) 5 a(i)(j) + w im S a(m)a) + J m S a(i)(m) , (4.5) 
-> hJJ - a k d k hj] - u> kl x l d k hJJ + uPhJj + u jk hj£ , (4.6) 

A label index was attached to the matter variables now. In (4.5) the transformation property of the local 
basis was used. Notice that the body-fixed basis in (4.5) stays unchanged. For (4.6) and (4.7) the transverse- 
traceless projection was commuted with the infinitesimal coordinate change. 

4.1.2 Symmetry Generators 

Now we try to construct the generators of infinitesimal rotations and translations, Pj and Jy,. These are 
of course nothing else than 3-dimensional total linear and angular momentum. With the help of these 
generators the transformation rule for an arbitrary phase space function A must read 

A -> A + \u ij {A, J ]t ) + a l {A,P t } . (4.8) 

It is sufficient to guarantee this transformation rule for all canonical variables. Comparing (4.8) with (4.4— 
4.7), using the standard Poisson brackets (3.48, 3.49) for each object as well as (3.62), one can indeed 
construct Pj and J^ . It is immediately clear that Pi and are a sum of matter and field parts, 

Pi = P 4 mattel ' + Pi dd , Jij = Jg"* + 4 dd • (4.9) 

In order to get Pi, one sets w l i = and a 1 is taken to be arbitrary. Then among the matter variables only z l a 
is transformed. Comparing (4.8) with (4.4) one obtains Sij — {z l a , Pj} — gp 2 - for each particle, and thus 

pmatter^^. (410) 

a 

Similarly, for the field part one gets —dkhJJ = {hJJ, Pf.} as well as —d^ 1 ^ = {7r yTT , Pfe}, which leads 
to 



1 

i&TT 



pMd = --±-\ d 3 a .^rPT^TT. (4 n) 



The derivation of Jij is analogous, with the result 

■or* = - + e . ( 4 - 12 > 

a a 

-_i-|d^2(^ TT ^-^ TT ^). 



(4.13) 



The ADM Hamiltonian Padm is by construction manifestly invariant under global rotations and trans- 
lations (at least in the considered gauges). Comparing Padm #adm with (4.8) one sees that total 
linear and angular momentum have vanishing Poisson brackets with the ADM Hamiltonian and are thus 
conserved. 

Yet another symmetry specific to objects with spin is given by constant rotations of the body-fixed 
frame, 

Af?^ -> A l * Kj) + oj^A^ , (4.14) 



4.1 Symmetries and Conserved Quantities 
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parametrized by a constant antisymmetric matrix ui a = -w„ 1J for each object. The corresponding gen- 
erators read = ^a^^-a^^S a ^)(i) and are also conserved quantities, ^^jui = const. A corollary 
of this is that 

jh aim J a°im = = COnSt ' ( 4 - 15 > 

As the ADM Hamiltonian Hadm is invariant under the transformations (4.4—4.7) and (4.14), these 
transformations are also a symmetry of the action (3.63). The corresponding conserved quantities Pi, Jy, 
and J?m can then be obtained by standard Noether arguments [90] and come out identical to above 
results. 

4.1.3 Global Poincare Invariance 

The global symmetry under the Euclidean group discussed in the last section is only a part of the bigger 
global symmetry under the Poincare group. Besides total linear and angular momentum, also the boost 
vector J l ° and the energy E — Hadm of the system generate a symmetry of the action and are conserved 
quantities for asymptotically fiat spacetimes. However, the infinitesimal transformations generated by 
Hadm and J l °, similar to (4.8), are in general highly nonlinear in the considered gauges and may not 
be written down immediately, as opposed to (4.4-4.7). Further, J i0 explicitly depends on time, see (2.12). 
Hadm and J i0 can be calculated by surface integrals at spatial infinity, see, e.g., [54, 91]. For the total 
energy E = Hadm this was already found in (2.45) and the boost vector J M is given by (2.12) with 

G i = — <j)d 2 s k [a; i (7 W)i - 7«, fe ) - 7fc + S lk j u ] . (4.16) 

Similarly, for 3-dimensional total linear and angular momentum it holds 



p i = -— ( t > d2g ^ lk , Jij = ~ZZf d 2 s fe (xV fe - x^ lk ) . (4.17) 



When these quantities are expressed in terms of canonical variables (after gauge fixing), they fulfill the 
Poincare algebra (2.10, 2.11). Notice that all Poisson brackets in (2.10, 2.11) involving Pi and just 
reflect the transformation property (4.8). Similar to the special relativistic case in section 2.1, one can define 
different total spins and centers for a gravitating system in asymptotically fiat spacetimes. In particular, a 
center and total spin of the system with standard Poisson brackets can be constructed (this was exploited 
recently in [92]). 

4.1.4 Symmetry Generators from Integral Formulas 

For simplicity we assume that 7y does not need a redefinition in order to receive a canonical meaning. 
However, this might be necessary at higher orders in spin. For the canonical field momentum we allow 
spin corrections by the ansatz 

7C ij = n ij + 16tt < j K , (4-18) 

a 

where n% contains the yet undetermined spin corrections. The gauge condition 7r M = with the subsequent 
decomposition (3.60) is assumed to hold. The surface integrals from the last section can be transformed 
into volume integrals using the Gauss theorem. With the decomposition (2.50) it follows 

E= \- { d 3 a;Af G i = [ d 3 xa; l A0. (4.19) 



167r J I6ir 

However, it is not possible to express E and G l in terms of the canonical variables without solving the 
nonlinear constraint equations for <j). Similarly one gets 



Pi = ~^j d 3 X7T lfc , fe , J ij = -l-[ dM^TT^fc 



x j 7T ik !k ) . (4.20) 
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Here one can exploit the momentum constraint Hi = 'Hf eld + %™ atter = to further evaluate Pj and Jy 
without needing to actually solve the constraints. Using (2.44, 2.50, 3.60) the momentum constraint can 
exactly be written as 

7r ifc ifc = -8ir{H™ mr + H~ T ) + ^ kJJ hy k l ~ (^ fcTT C),j - A (V k h^ + B ij tj , (4.21) 



with the definitions 

^matter = £ [^.^ _ 2 ( 7lfc n k J5 a ) , 



(4.22) 



B« = [l - (1 + |0) 4 ] (#« + tt« tt ) + . + ^ - - ^,fe^ T , (4.23) 

and the alternative vector potential 

V* = (Su - -didjA-A , (4.24) 



4 

for which it holds 

* ii = V i lj + V i , i -~6 ij V k , k . (4.25) 

One can calculate y,f Mer using (3.41). Notice that B ij = B 3i and B u = 0. Further the last two terms in 
(4.21) do not contribute to (4.20). Obviously, (4.20) are a sum of matter and field parts, (4.9). The field 
parts are identical to (4. 1 1) and (4. 13). However, the matter parts now read 

pmatter = f ^ x ^matter + ^matter) ^ (4 26) 

jmatter = f g x (^matter + ^jmatter _ ^matter _ ^mnatter) _ (4 2?) 

For consistency, these must turn into (4. 10) and (4. 12) when expressed in terms of canonical variables. 

4.2 Construction of Canonical Variables 

In this section we will formulate the important consistency conditions and apply them order-by-order in 
the post-Newtonian approximation to find canonical variables. 

4.2.1 Consistency Conditions 

In section 2.2.3 it was seen that the spin length S a given by 2S% = S a [ lu S% u is a conserved quantity 
in the covariant spin supplementary condition. This can also be derived from the action (3.30) using 
the symmetry under constant 4-dimensional Lorentz transformations of the body-fixed frame, see also 
[69]. This conserved quantity must be identical to the one in (4.15), as both were derived from the same 
symmetry (though only the 3-dimensional rotation part is relevant after the supplementary conditions were 
eliminated). Thus it must hold 

SafiuS^ — Sa(i)(j)Sa(i)(j) , (4.28) 

providing a relation between covariant spin S aflll and canonical spin S a ({^/j^ . This is one important consis- 
tency condition we will impose. 

Further, one can calculate "H™ atter and thus (4.26, 4.27) in terms of (noncanonical) variables in the 
covariant supplementary condition with the help of (3.41) and (2.30), 



-^matter ^ ^ 



(p ai + K t] nSl)5 a + ( \^ k S aik 5 a + yk PadnSak) g \ 



(4.29) 
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see (3.40). Then (4.26, 4.27) must coincide with (4.10, 4.12), leading to conditions on the transformation 
between canonical variables and variables in the covariant supplementary condition. We write this as a 
condition on %™ attel in the form 



(4.30) 



where the symmetric part of s l J is not constrained, but it has to hold 



S, 



(4.31) 



This condition on "H™ attel is the most general one 13 that guarantees that (4.26, 4.27) coincide with (4.10, 
4.12). 

Above conditions are sufficient for the post-Newtonian order considered here. Another condition that 
could be useful at even higher orders (in particular also higher orders in spin) would be the fulfillment of 
the Poincare algebra. However, all Poisson brackets in (2.10, 2.11) involving Pi and are fulfilled by 
construction due to the transformation property (4.8) if above conditions hold, thus giving nothing new. In 
[14] it was considered whether the construction of the constraint algebra (2.72-2.74), which is related to 
diffeomorphism invariance and thus more fundamental than global Poincare invariance, could be used to 
construct canonical variables. However, this approach seems to be unmanageable. 



4.2.2 Canonical Variables 

First we evaluate the condition on the spin length given by (4.28). We will first construct a specific trans- 
formation between S a ij and S a ^\fj\ and then discuss its uniqueness. Inspired by the flat space case (2.20), 
we first apply the transformation 



Pai n S a 



PajTlS 

it 



np a 



to the conserved quantity S atll/ S£ u 



, nS a i 
- np a 

l kl l l] S akl S ai3 



Pakl 



kj c 

•3(1 



(4.32) 



2 1 ll nS ai nS aj , with the result S aflv S^ v 



ry«yj SakiSaij- With the help of an arbitrary triad eu\j this can be written in a local basis as S a ^ v S^ v = 
S a a)U)S a u\fj\, so we have found a possible transformation allowed by (4.28). The ambiguities that are 
left can best be discussed in terms of the spin vector S a ay As we are still considering the linear order in 
spin, any further transformation of S a {i) must be linear in spin and must leave the expression S a u)S a <i) 
invariant (notice S a (i)(j)S a (i)(j) = 2S a ^S a ^). Therefore only a rotation of the spin vector as a further 
transformation is possible, which can be absorbed into the yet arbitrary triad e^j. 
A comparison of (4.29) with (4.30) leads to 



Pa 



Pa 



K i3 nSi + 7rf 7j - M 



-Sakj 



Pa(k' 



iS, 



a.3) 



np a 



(4.33) 



without any ambiguity. Now (4.29) is of the form (4.30), so (4.31) is the only condition that is left. In order 
to evaluate (4.31) we first need to read off sjf. For the redefinition of the position variable we use 



np a 



'An ) 



(4.34) 



where z\ a is a yet unknown correction to the flat space case (2.20). Comparing (4.29) expressed in terms 
of the new variables with (4.30) leads to 



in ik A ik In ^PalPaii^ ak)p „ A j 

s l i = y k S atk + 7 3 7 P s ( " 2 P<^A a , 

npa{m a -np a ) 



(4.35) 



13 In the pole-dipole approximation at most one partial derivative can appear in %™ altel . Further it was assumed that the variables 
from different objects do not mix (e.g., as P1&2) at this stage. 
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with the definition np a = — -J m 2 + J^paiPaj ■ The only ambiguities in the transition to canonical vari- 
ables are now given by ir l J , z l Aa , and the triad eu\j. We try to fix these ambiguities by considering (4.31) 
with (4.35) order-by-order in the post-Newtonian approximation, which is introduced in the next section. 

From the action approach we know that the ambiguity of eusj should just be a gauge freedom. Thus 
different choices for eu\j should be canonically equivalent. Indeed, it was shown in [84] that a spin rotation 
is just a canonical transformation at linear order in spin. However, a canonical transformation may change 
all variables, but p a i as well as hJJ can not be changed any more. Thus the canonical representation was 
already partly fixed and we must therefore still keep e^j as general as allowed by the restriction on the 
triad gauge made in section 4.1.1. 

4.2.3 Post-Newtonian Approximation 

The idea behind the post-Newtonian approximation is that for slowly moving bodies and weak gravitational 
forces the Newtonian physics is recovered as a first approximation. For two objects this means that their 
relative velocity v shall be small compared to the speed of light c. In Newtonian physics the time average 
of kinetic and potential energy is of the same order if the virial theorem applies, which is the case for bound 
systems. Then one has 



where M is the total mass of the system and r the typical distance of the objects. An expansion in the 
dimensionless quantities (4.36) obviously is also an expansion in c~ 2 . We will therefore think of the post- 
Newtonian expansion as an expansion in c -2 . However, this is a rather formal point of view as it depends 
on the choice of units whether c~ 2 is actually a small number (e.g., in our units it is equal to one). As seen 
later, there may be half post-Newtonian orders corresponding to c _1 . 

As post-Newtonian orders are formally counted in terms of the velocity of light c originally present 
in the equations, i.e., before setting c = 1 = G, one should introduce G and c back into all expressions. 
However, this would undo the advantages achieved by setting c = 1 = G. Instead, we will assign an order 
in powers of c _1 directly to our variables. When setting c = 1 = G only one unit is needed, which we 
choose to be the unit of spatial distances, e.g., meters. Then the values of all masses m a must be given in 
meters, which is obtained by multiplying their values in kilograms by G/c 2 . Therefore we just count the 
masses to be of the order c~ 2 , as this is the power of the speed of light that would be introduced into the 
expressions if we restore the original units. Similar arguments apply to the other matter variables and we 
have the counting rules 



Notice that an energy receives a counting of c~ 4 , which gives the absolute order of the Newtonian Hamil- 
tonian (being an energy) within these counting rules. However, one obtains different counting rules for 
the matter variables if one uses kilograms instead of meters to replace all units when setting c = 1 = G. 
This convention is also often used and leads to different absolute orders in c _1 , e.g., a mass now receives a 
counting of c° and the Newtonian Hamiltonian is at the absolute order cT 2 . But relative orders are always 
the same, so only a counting relative to the Newtonian order (or to the leading order if the Newtonian order 
vanishes) finally makes sense when using such counting rules. The correct absolute Newtonian order is c°, 
as it must prevail when c _1 — s- 0. 

The formal counting may nicely be applied to more complicated situations, e.g., when spins are present. 
For dimensional reasons only we are thus counting the spins of the order c~ 3 . This has some computa- 
tional advantages, e.g., similarities to calculations for nonspinning objects are more manifest, see section 
6.2.2. Here post-Newtonian orders should always be understood in the formal sense if not otherwise stated. 
However, the spin of a (Kerr) black hole is given by Gm 2 a/c, where m is the mass of the black hole and 
a = ... 1 is the dimensionless Kerr parameter. The maximal spin of an object is defined as Gm 2 /c 
(which is the maximal spin of a black hole, a = 1), and additionally has to be counted as c _1 . If the spins 
are maximal, one therefore has to add half a post-Newtonian order relative to the formal counting for each 
spin variable appearing in a specific expression. 






(4.37) 



4.2 Construction of Canonical Variables 



31 



If the spins are not maximal, one has to be careful when classifying spin effects into post-Newtonian 
orders. For example, if the spin is of the maximal one and the orbital velocity is of the speed of 
light, then each spin variable corresponds to one extra order in v/c relative to the maximal spin case, or half 
a post-Newtonian order. At a later time during the inspiral the spin length has not changed much 14 , however, 
the orbital velocity might have increased, e.g., to A of the speed of light. Then each spin variable even 
corresponds to two additional orders in the velocity or one post-Newtonian order relative to the maximal 
spin case. To conclude, while the spin length does essentially stay constant during the inspiral, the orbital 
velocity will increase and one expects that spin effects slightly shift to higher post-Newtonian orders during 
inspiral. Therefore, assigning a post-Newtonian order to spin contributions in the Hamiltonian seems to 
make no sense in general, except for maximal spins or within the formal counting. However, this discussion 
is only superficial, the relevance of spin effects also crucially depends on the orientation of the spins and 
the mass ratio of the objects. Due to these problems we will prefer to classify spin effects by leading order, 
next-to-leading order, etc. when possible. 

Counting rules for other quantities may be derived from (4.37). For example, cj> results from solving 
the constraints, and one may easily see that its leading order must be identical to the leading order of the 
matter source of the Hamilton constraint "H matter , which is cr 2 (this will become obvious in section 6.1.2). 
Similarly one gets counting rules for the other field variables by considering the matter source of the field 
equations. Without going into detail, we state here that 

4> = O (c~ 2 ) , hjj = O (c- 4 ) , jf y = O (c- 3 ) , n ijTT = O (c~ 5 ) , (4.38) 

are the correct counting rules for the fields. In general the fields include different post-Newtonian orders, 
(4.38) only gives the leading orders. The Taylor expansion of the fields in terms of c~ 4 is written as, e.g., 

- 0(2) + 0(4) + 0(6) + O (c- 7 ) , (4.39) 

where a number in round brackets denotes the absolute order in c _1 within the counting given by (4.37) 
(this should not be confused with indices in the local basis). The vanishing of the odd orders <p^ and <j>/g\ 
is explained by the vanishing of the corresponding orders in the source terms. 



4.2.4 Final Fixation of the Canonical Variables 

First we try to find a way to parametrize the ambiguity in the triad when the induced metric is kept fixed. 
If one considers the perturbative expansion of e'^'e 1 '''' = ^ under the assumption that the leading order 
is given by e^? — 8ik, then one sees that the symmetric part of e 1 ^ is uniquely fixed at each order, while 

the antisymmetric part e lJ = |(e*^ — e J w) is arbitrary. Therefore e iJ parametrizes the rotational degrees 
of freedom left in the definition of the local basis and thus the ambiguity of the canonical spin variable. In 
particular, the leading post-Newtonian orders read 

4f - - jMw . 4f - % - 5«§)$ " jMw + !M(2) - \h% . (4.40) 

Notice that e 4 - 7 is needed only on the worldlines. In the following we use the abbreviation e l J = e lJ 

Next we make an ansatz for ir l J , z Aa , and e z J at each post-Newtonian order. For this purpose it is 
important that ir l J has the dimension length squared, z l Aa the dimension length, and e l J is dimensionless. 
Further, 7rJ? and z Aa must be linear in spin, while e l J must be independent of the spins. The fields hjj and 
Ti-y'TT are a j wa y S taken at the position z l a in such an ansatz and z l a should not appear directly. Of course 
one also has to take into account that tt' 1 J must be symmetric and antisymmetric. Considering possible 
ansatze under these restrictions we infer that the leading orders are at least 7T* J = O (c~ 5 ) , z Aa = O (c~ 2 ) , 

and e'J — O (c~ 6 ). From (4.35) the first orders of Sa^ then follow as 

S L(3) = Sa(i)(j) > s i(5) = PajZ\ a {2) ~ Pai z Aa(2) ' S a(7) = ^aj 4a(4) _ Pai z A a (i) ' ( 4 -41) 



In the approximation considered here the spin length is even exactly constant. 
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Evaluating (4.31) one concludes that zXa(2) = an< ^ z Aa{4) ~ 0- Thus we have anticipated the correct 
redefinition of the position (4.34) to this order. 
For s^/L one has 

o(9) 

S a(9) = g a fc (6)^a(fc)0) " Pa<4a(6) + ^2?^ {p al S a(l)(k) + p a lS a{l)(k) ) j) , (4.42) 



where (i <H> j) denotes an exchange of the indices i and j. The most general solution of (4.31) under above 
restrictions is 

1 — C 

n a{6) = -^r(P^kSa(k)U) + PajPakS a{k){l) ) , (4.43) 

G C 

g a(6) = -^PakiPaihJl - p a jhj£) , zXa(6) = J^2^« ( S a(k)(i) hjk + S a{k){j )hJl) , (4.44) 

at this order and now depends on an arbitrary constant C. 

However, we can remove the ambiguity C by a canonical transformation with an infinitesimal generator 



C 

9 = -^—2paiPakS a (k)(j) I d A xh\-8 a . (4.45) 



An arbitrary phase space function A then transforms as A — >• A + {A, g) to the required order. Applied to 
the canonical variables one obtains 

hjf^hjj, -j- ^ TT - 5™ 3 X ^P ofc p am 5 a(m)(0 5 o , (4.46) 

a m<1 

c 

C ~> - ^Aa(6) 5 Pai^Pai- J^PalPajSaVXtyhuj. (4.48) 



This indeed removes all terms depending on C from the source expressions % mattel and %™ atter at the 
considered order. We can therefore choose C = 0, which leads to agreement with the action approach. 
The triad then is in the spatial symmetric gauge e iJ = and all variable transformations are the same as in 
the action approach at the considered post-Newtonian order. In particular, using (3.55) and (3.58) in (3.45) 
leads to (4.33). Further, in the action approach we found that z l Aa = and 

^ = \A^+B%Af\ (4.49) 
or more explicitly using (3.44) 



_/./ _ „.ik„,jl m aPa(k n Sal) ^ ^ffj^km^ln 

2np a (m a - np a ) 2 " 



7*V r ^ + t: B k h km -y ln S amn . (4.50) 



Using BV, = O (c~ 4 ), cf. Eq. (3.37) in [15], the post-Newtonian expansion of (4.50) agrees with the 
findings in this section. The check of the action approach given here is valid to the formal 3.5 post- 
Newtonian order. 



5 Higher Orders in Spin 

Higher orders in spin require higher multipole moments, e.g., a black hole has a nonzero quadrupole at the 
quadratic level in spin [37]. We will constrain to quadrupole and quadratic order in spin in this section. 
Besides spin-induced quadrupole deformations discussed here, also tidal deformations induced through the 
gravitational field of other objects have been treated in the literature, see, e.g., [93, 94]. 



5. 1 Quadrupole Approximation 
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5.1 Quadrupole Approximation 

The extension of the pole-dipole approximation to higher multipoles was already essentially completed 
some time ago [41, 95-97], see also [98, 99], most notably by Dixon. It should be stressed that Dixon's 
method incorporates Mathisson's pioneering ideas [100]. 



5.1.1 Quadrupole Approximation from Tulczyjew's Method 

A more direct application of Mathisson's ideas to the quadrupole order was given in [16] with the help of 
W. M. Tulczyjew's method [31], see also [101]. There the quadrupole moment t^ va ^ was kept in (2.26) 
when evaluating (2.27). In addition to p p and S^ u now various quadrupole moments appear. It is suitable 
to introduce a reduced quadrupole moment J pua P with symmetries 

jvp/3a _ j[vp][f3a] _ jfiuvp ju[pf)a] _q ^ jvpfiu _^ jvfiap ^_ jvapfi _ g (5 1) 

Thus JpP av has the same (algebraic) symmetries as the Riemann tensor. This quadrupole moment is able to 
incorporate all quadrupole contributions from t pua ^ that remain after (2.27) was evaluated. The equations 
of motion then take on the simple form 

— = 2p^ + -R (A l ["J"]'** , ^ = --R?\ u p S? a - -R { % „ r pf3a , (5.2) 
3 al3p dr 2 pp<3a 6 vp P a W p 

and agree with Dixon [97]. The reduced moment J pva P is also optimal to give a simplified expression for 
the stress-energy tensor, reading 



~gT pv = / dr 



(5.3) 



In this form the stress-energy tensor was first given in [16]. This stress-energy tensor, together with the 
ansatz for J pua at the quadratic level in spin given below, can be applied to the derivation of the next-to- 
leading order radiation field, see [102] for the spin-orbit case (the leading order is given in [103]). Besides 
this formula for the stress-energy tensor, a further interesting result in [16] is the relation between the 
£/«/... momen t s anc j Dixon's reduced moments p p , S pv , and J pva P . This relation could be used to study 
alternatives to Dixon's integral formulas for the multipole moments or to discuss the relation between 
moments belonging to different representative worldlines (for the latter see section VII. C in [16]). 
The spin supplementary condition S pv ' f v = is preserved in time if 



which should give a relation between p p and u p . This extends (2.31) to the quadrupole approximation. 
The extension of (2.35) reads 

p" = -u^u p ^pu„ + f u^/j^*" . (5.5) 



5.1.2 Decomposition of the Quadrupole 

In order to parametrize the quadrupole deformation due to spin we try to find the most general covariant 
ansatz for J p - va P quadratic in the spin tensor that is relevant for the post-Newtonian order in question. It is 
suitable to consider the orthogonal decomposition of J pva P with respect to the vector to which the spin 
is orthogonal, S pu f v = 0. This decomposition reads 

jvpPa = Qvpfia _ 1 (fWQP]P<* + f[<*Q0]P») fW QP]\P f«] ; (5.6) 

where Q^ Q , Q p/3a , and Q pf} are called stress, flow, and mass quadrupole here and are orthogonal to 
in each index, see also [104]. They further have the symmetries 

Q V pfj a = Q[up][pa] = Qfiavp ( Qu\pPa] =Q ^ Qvptia + Qvj5 a p + Qvapfi = q ^ ^ ?) 
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Q P p a = Q P [0a] ^ QlpfSa] _ g ^ Q P p a + Q0ap + Qap/3 = g ^ gp/3 = g(p0) _ (5 g) 

In a local basis with / M giving the time direction these moments only have spatial components (due to 
the orthogonality of these moments to f p ). One may therefore decompose these moments further in the 
local basis into parts transforming under irreducible representations of the 3 -dimensional rotation group 
SO(3). For the mass quadrupole this SO(3)-decomposition reads, written in the coordinate frame, 

Qfj.u = Q pu ^ T^PpvQ^ p i (5-9) 

with the orthogonal projector P pv = g pv — yypf >J 'f v (notice P pv ' g py = 3). Here Q S JJ is symmetric 
and trace-free (STF) in the local frame. In the coordinate frame the trace-free property reads Q^Jg^ = 
and of course it holds Q^JJ" ' f v = 0. Obviously Q p,v has six independent components, five contained in the 
symmetric trace-free part and one in the scalar part Q p p . The same holds for Q p " a P with a more compli- 
cated decomposition into symmetric trace-free and scalar parts, whereas Q pva has even eight independent 
components corresponding to a symmetric trace-free and a vector part [104]. 

5.1.3 Ansatz for the Mass Quadrupole 

We will now constrain to = p p and to the Newtonian limit, the latter to identify the dominant contribu- 
tions. Then only the mass multipoles are important for the dynamics. Though flow and stress multipoles 
do in general not vanish in the Newtonian limit [96], they give no contribution to the gravitational field and 
can be neglected. Therefore the decomposition of the quadrupole moment (5.6) just reads 

jupf} a = _ 3 [uQp][P p a] (5 1Q) 

m p 

where the dynamical mass defined by PpP p = —m p is now denoted as m p . Also the trace part of the 
mass quadrupole gives no contribution to the gravitational field outside the body and one can thus assume 
Q p p = 0. The mass quadrupole induced by spin is then given by the ansatz 

Qfj,u = Q fn, — — ( Sfj_pS u p — 7;PpijS p Sp a J , (5.11) 
nip \ o J 

and is parametrized only by Cq in the Newtonian limit and quadratic level in spin, see also [105]. For 
black holes one has Cq = 1 [37] while for neutron star models Cq depends on the equation of state [106]. 

Though the ansatz for the quadrupole (5.10, 5.1 1) was given in the Newtonian limit only, it was written 
in a manifestly covariant way and we will now consider its implications in full general relativity. However, 
we stay at the quadratic level in spin. It will be shown in the next section by relying on investigations in 
[107] that this ansatz indeed holds to next-to-leading order in the post-Newtonian approximation. It is easy 
to see from (5.2) that the spin length S given by 2S 2 — S^S^ is conserved for our quadrupole ansatz and 
spin supplementary condition S pv pp = 0. But the mass m p is not conserved. However, the new mass-like 
parameter m defined by 

™ = m P -lRi% a r p0a , (5.12) 
is conserved for our ansatz quadratic in spin. Finally (5.4) can be written as 

2m Kp ^ b b U + 2 
and gives a relation between p p and u 



5.2 Action Approach 

It is shown in this section that allowing nonminimal couplings in the action approach from section 3.1.3 
corresponds to certain higher multipole corrections, see also [87]. The couplings in the action needed for 
spin-induced quadrupole deformations at next-to-leading order in the post-Newtonian approximation can 
be found in [107]. 



5.2 Action Approach 
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5.2.1 Nonminimal Couplings 

We now generalize the ansatz (3.16) for the action to nonminimal couplings. More precisely, the La- 
grangian is allowed to contain the Riemann curvature tensor, 

All _ I A^T.J*,V OM" „&v(„P\ n („P\ E>( 4 ) 



W M [ei^z^K A1 ] = J dTL M K,^,r(/),v(^a<^M)- (5-14) 

The Euler-Lagrange equations of this action follow as in section 3.1.3 in a straightforward way. It is easy 
to see that (3.19) stays unchanged, 



D 

d7 



dL 



M 



^Lm^p ^^m^p (515) 



However, the important relation (3.17) now reads 

U ~ du° + dfi- + dg«» g dgp v 9av d R^ s " 

Using this identity and the definitions 

the Euler-Lagrange equations for the matter variables turn into (5.2), whereas for the field variables one ob- 
tains the Einstein equations with the stress-energy tensor (5.3). Higher multipoles are covered by including 
symmetrized covariant derivatives of the curvature tensor in the Lagrangian [87]. 

An action invariant under general coordinate transformations always leads to a stress-energy tensor 
fulfilling (2.27). As well known, this can be shown from the Noether identity [90] following from general 
covariance, see also Eq. (18.23) in [108]. However, this does not mean that an action approach as envisaged 
here always leads to the most general equations of motion allowed by (2.27). Dixon's derivation essentially 
only evaluated (2.27) and thus covers a much more general situation, though the Euler-Lagrange equations 
obtained in this section are identical to the equations of motion found by Dixon (but it is not a priori clear 
that this will be the case). In particular, the multipole moments will always be implicitly fixed by the other 
variables on which the Lagrangian was chosen to depend on, cf. (5.17). This means that the quadrupole 
is not a dynamical variable within our action approach and our action does not cover, e.g., quadrupole 
oscillation modes or tidal resonances, see, e.g., [109]. Notice that (2.27) puts no constraints on equations 
of motion related to dynamical quadrupole degrees of freedom. However, a dynamical quadrupole requires 
further dynamical variables in the action principle. For a good effective description of extended objects via 
an action one thus needs some intuition on the relevant degrees of freedom that should enter into an ansatz 
for the effective action. The action (5.14) includes translational and rotational degrees of freedom of the 
object, which is expected to be a good choice for the inspiral phase. 

An advantage of the action approach is that it is easier to find conserved quantities or constant pa- 
rameters. In particular, the conservation of the spin length immediately follows from the symmetry under 
Lorentz transformations of the body-fixed frame. Further all parameters in the action, e.g., a mass-like 
parameter, are constant simply by assumption. It is much more difficult to find such constant quantities if 
one only considers Dixon's results together with a specific ansatz for the quadrupole moment, see, e.g., Eq. 
(5.12) or the discussion in reference [16]. 

It is important that one may eliminate the Ricci tensor (and scalar) from the matter Lagrangian Lm 
by a suitable redefinition of the metric [110]. This was already found in [111] within a slightly different 
situation and is based on the observation that in a perturbative context the use of lower order equations of 
motion in the perturbation part of the action corresponds to a redefinition of variables, see [112]. If we take 
some additions to the point-mass Lagrangian as a perturbation, then we may eliminate the Ricci tensor by 
using the Einstein field equations with the point-mass stress-energy tensor as a source, corresponding to an 
irrelevant redefinition of the metric. However, the point-mass stress-energy tensor then gives rise to singular 
self- interactions in the matter Lagrangian Lm, which are formally neglected [111]. The conclusion is that 
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one may use the vacuum field equations Rj$ = in the matter Lagrangian Lm- (This will also be used in 
a slightly modified way in section 5.2.3.) The matter Lagrangian can therefore be restricted to depend on 
the completely trace-free Weyl tensor C < ^J aj3 , 

c tlv(s = r( ^Lp + 9 a [vR { ^ - g^Rf^ + -^9p,[ugi3] a R (i) , (5.18) 

(4) 

instead of R\J, a fi' ^his wou ld gi ye r i se to corresponding modified multipole moments defined analogous 
to (5.17). Further, the Weyl tensor can be split into electric E^J and magnetic B$ parts, 

B (4)_ (7 (4) a a d(4) _ 1 (4) c (4) a B p a ,c m 

with tj^pp the 4-dimensional Levi-Civita symbol, leading to definitions for corresponding electric and 
magnetic multipoles as partial derivatives of Lm- It could be interesting to consider the impact on the 
equations of motion and the stress-energy tensor from letting Lm depend on C^J a/3 or E$ and B$ 
instead of R^^p- However, this will not be necessary here. 

5.2.2 Legendre Transforms and Supplementary Conditions 

Up to now the Lagrangian Lm is completely arbitrary and the equations of motion fully agree with Dixon 
at the quadrupole level. The question is which supplementary conditions (3.10) and (3.12) belong to Lm, 
or how Lm must be chosen to fit with specific supplementary conditions. We only require here that (3.10) 
and (3.12) are preserved in time, which leads to 

+ A AI (r,„ - j^j ^ = , S%W»f v + = , (5.20) 

where (5.15) with (5.17) was used. Both conditions are fulfilled if we have 

Q^fu + P^— = 0- (5.21) 
dr 

In this sense the conditions (3.10) and (3.12) belong together (however, there may be exceptions). This is 
the condition our action shall fulfill here. Notice that (5.4) only guarantees that (3.12) is preserved in time 
so that the second relation in (5.20) holds, but this does not imply (5.21). However, comparing (5.17) with 
(5.4) can still be useful. Further, at the quadratic level in spin we need to fulfill (5.21) only to linear order 
in spin. 

It is suitable to define a new function Rm(u^, S^ v , g^ v , g^ v , R^„ a p) via Legendre transformation, 
R M = L M - ^S^W. It holds 

An ansatz for Rm then has to fulfill the condition (5.21) with (5.22) inserted. This gives a partial differential 
equation for Rm- It holds Rm = P^u^, which is a consequence of the reparametrization invariance of the 
matter action. Notice that Rm is similar to the Routhian used in [107, 113, 1 14]. For the Routhian the Ricci 
rotation part in the term hS^vW, cf. Eq. (3. 15), is not subtracted from the Lagrangian Lm- Therefore the 
Routhian is not a covariant function, whereas Rm introduced here is covariant. 

Due to reparametrization invariance a full Legendre transformation in it M and leads to a vanishing 
result. However, as in section 2.3.2 we define a function Hmt which contains the mass-shell constraint, 
suitably generalized to the quadratic-in-spin level, together with a Lagrange multiplier A. It holds 

U » = ?EML ^ = 2 d^r j^ = 6 9Hm^ (5 23) 

dp, dS,„ BR W ap 
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It is also possible to give an ansatz for the mass-shell constraint and thus for Hmt directly. This ansatz 
must be chosen such that the condition (5.21) with (5.23) inserted is fulfilled. Further, one may simplify the 
quadratic-in-spin corrections to H m t by using the leading order constraint p^p p = — m 2 , corresponding 
to a redefinition of the Lagrange multiplier [112]. 



5.2.3 Leading Order 

The coupling terms found in [107] adapted to our notation and conventions read 

R M = . 1 = fm«X - ±R^ S^S^u»u p + 9l^Ef}S%sA . (5.24) 
y/ —u a u a \ 2m p 1 2m p / 

It was found in [107] that these coupling terms are the most general ones at quadratic level in spin sufficient 
for the next-to-leading order in the post-Newtonian approximation. Besides these terms corresponding to 
quadrupole deformation due to spin, one could also treat tidal deformations, see, e.g., [93], using nonmin- 
imal couplings in the action given in [94, 115]. With the equivalence of Riemann and Weyl tensors within 
the matter action, see section 5.2.1, we can write Rm as 

Rm = 7=7^ { mU ^ ~ ^ R ^S^S^u v u p - ^i^Q^u'v) , (5.25) 
where Q^ v is given by (5.1 1) and Ces 2 = Cq. If we set f ^ = p M and thus 

S»» Pv = S^—^=0, (5.26) 

we obviously reproduce (5.10) and (5.13) within the gauge u a u a = —1 by plugging (5.25) into (5.22). 
Further (5.21) is fulfilled 15 to the considered order in spin by using (5.2) and (5.22), i.e., 

n - g di? M dR M _ 1 ptla (4) p „ fia p/Mr (4) dR M 
* " L upf3a 

The last term is of higher order here as this condition must be fulfilled to linear order in spin only. Notice 
that m depends on spin according to m — too + jjS a pS a ^ , see also (3.28). (Otherwise the Legendre 
transformation between Lm and Rm would not be possible.) 
An equivalent description in terms of Hmt reads 

Hm t = A ^to 2 + Vv3 r + -L R ^ a0 S^S^p» Pp - ^^S"^) , (5.28) 

with the action still given by (3.30). The derivation of a canonical formalism now follows along the same 
lines as in section 3. First the matter constraints are solved. The only difference to the linear-in-spin case 
arises in the mass-shell constraint, which follows from the variation of A. The solution of this constraint 
reads 

up = n»p, = -J m 2 +1 H PiPj + ^ <^S%S P Vy • (5.29) 

The last term was not yet split into time and space parts. This splitting leads to quite many terms, so one 
should restrict to some post-Newtonian order. Though all formulas are sufficient for the next-to-leading 
order, we will for simplicity only treat the leading order in this section. Then we have 

np = - v / m 2 +7 o- piPj ._ _g_ 7 fc' 7 -y»5 ifc 5 j , (x mn , + N, mn ) . (5.30) 

The only contribution to the action quadratic in spin then arises from the term Nnp in the matter La- 
grangian, see (2.68). Problematic is the partial time derivative of the extrinsic curvature. In consideration 



15 One could also consider the most general ansatz for Rm and ask for which choice of the condition (5.21) is fulfilled. This 
would allow one to study the impact of the supplementary conditions on the dynamics. 
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of the definition 2iVify = — 7ij,o + ^^(i-.j) we see that the ^mn.0" term produces time-derivatives of 
lapse and shift, as well as a double time-derivative of 7^ . This does not fit well to the derivation of the 
canonical formalism as given in section 3. In order to overcome these problems, we eliminate K mn< Q with 
the help of the vacuum field equations, cf. the discussion in section 5.2.1. Finally one ends up with just a 
quadratic-in-spin correction "H™" 61 to the source of the Hamilton constraint % matter of the form 



This source term is quite unusual in the sense that it is not a specific projection of the stress-energy tensor 
(5.3, 5.10, 5.11), i.e., % matter ^ \F(T llv vP"n 1 ' . This is due to the implicit redefinition of variables performed 
by using the vacuum field equations in the matter action. However, the leading order Hamiltonian resulting 
from this source term is identical to the well-known one obtained in section 6.2.2. The variable redefinitions 
from section 3.2.2 are still correct at the leading order. (There are no additional terms that need to be 
cancelled in the action and all quadratic spin contributions from the redefinitions in section 3.2.2 are of 
higher order.) 

At the next-to-leading order the calculation gets much more involved. In particular there are more time 
derivatives of the extrinsic curvature that must be eliminated and the variable redefinitions from section 
3.2.2 need corrections quadratic in spin. It is also relevant whether the field variables in the variable 
transformations are taken at the new or at the old particle position. Further corrections to the canonical 
field momentum seem to be necessary, too. We will therefore study an alternative derivation oriented at the 
symmetry generator approach from section 4 in the following. 

5.3 Symmetry Generator Approach 

We now sketch the derivation of the canonical formalism at quadratic level in spin via the approach from 
section 4. However, essentially only the calculation of the source terms of the constraints as certain pro- 
jections of the stress-energy tensor is used here, the determination of canonical variables by looking at the 
symmetry generators will only be touched lightly. 

5.3.1 Leading Order 



First we calculate the source of the field constraints as certain projections of the stress-energy tensor (5.3, 
5.10, 5.11), e.g., -H matter = ^fyT^n^n". To leading order we have 



and no corrections appear in "H™ attel . The variable redefinitions found at the linear order in spin are therefore 
sufficient here, as they are followed from % mat,er in the symmetry generator approach. For Cq a = 1 this 
source term is in agreement with the source of the Kerr metric in approximate ADM coordinates found in 
[116]. It further gives the correct Hamiltonian, see section 6.2.2. 

Obviously the derivation of the leading order in this section is much simpler than the one via the action 
approach. But this does not need to be true at the next-to-leading order. The problem is that it is not 
guaranteed that the variable redefinitions can be uniquely fixed by just the conditions (4.30) and (4.31). 
The action approach is much more systematic and should therefore be preferred at the next-to-leading 
order. However, in the next section a shortcut to the next-to -leading order Hamiltonian is described, which 
combines the approach of the present section with the Poincare algebra approach in [117]. 

5.3.2 Next-to-Leading Order Static Source Terms 

In [1 17] Hergt and Schafer constructed the part of the next-to-leading order Hamiltonian that depends on pi 
(i.e., the nonstatic part) from an ansatz for this Hamiltonian (together with a suitable ansatz for the source 
of the constraints). The coefficients in this ansatz could be uniquely fixed up to a canonical transformation 
by considering the Poincare algebra (2.10, 2.11). The degrees of freedom corresponding to the ambiguity 



matter 



2m 



(5.31) 




(5.32) 
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in the canonical representation are given by the coefficients that enter via an ansatz for the center of mass 
vector Gi. However, the static (i.e., f>i = 0) part of the Hamiltonian is left completely undetermined by the 
Poincare algebra approach in [117]. 

In order to get the complete next-to-leading order Hamiltonian only the static part of the Hamiltonian is 
missing, as well as the corresponding center of mass vector Gi. The latter is needed to consistently fix the 
canonical representation of the nonstatic part of the Hamiltonian given in [117]. Fortunately the center of 
mass vector does not depend on at the considered order. Therefore both the static part of the Hamiltonian 
and the center of mass vector are determined if we only know the static part of the source of the constraints. 
For pi — we get from the stress-energy tensor (5.3, 5. 10, 5. 1 1) 

KC:=o=Y,(h k ^ lJQ ^ s °) > (533) 



■kl 



but no further contributions to %™ attel arise. Though there is no difference to (5.32), this source term is now 
valid to next-to-leading order for the case pi = 0. 

However, we need the source terms for the case pi = and not for pi = 0. Also position and spin 
variables are not yet the canonical ones and we must discuss whether the variable redefinitions will have an 
impact on the source terms in the static case. As there are no contributions to %™ ttel ' at quadratic level in 
spin for pi = 0, no further static contributions to the redefinition of spin and momentum variables can arise 
from the conditions (4.30) and (4.31). Though static contributions to z Aa ^ could be necessary, they can 

be removed by a canonical transformation with generator —piz\ a uy (Notice that in the case pi = this 
transformation only changes the position variable.) Finally only the redefinitions found at the linear order 
in spin are relevant and only (4.33) gives contributions in the static case. The result for the static source 
finally reads 



a L \ / -,kl a 



:U (5.34) 



where 

Qaij = — — ( SaikSajl ~ xTy'T^T™" SakmSaln ) • (5.35) 
777 \ O / 

Equation (5.34) was found for the black hole case CQ a = 1 in [17] from a 3-dimensional covariant ansatz 
for Hp a * t !!f containing four coefficients. Two of these coefficients were fixed by matching to the Kerr 
metric, but the other two gave no contribution to the Hamiltonian or to the center of mass vector. One of 
the latter two coefficients would also arise here if we would have kept the trace part of the mass quadrupole, 
Q p p - The ansatz in [17] was generalized to arbitrary Cq a in [19]. 

The derivation given in this section is quite involved and it would thus be desirable to give a more coher- 
ent one with the help of the action approach in the future. This would also facilitate further investigations 
of quadrupole or higher multipole effects with the help of canonical methods. 



6 Results for Hamiltonians 

In this section the obtained canonical formalism is applied to calculations within the post-Newtonian ap- 
proximation. In particular, the next-to-leading order spin corrections to the conservative Hamiltonian are 
derived. The Hamiltonians are checked with the help of the global Poincare algebra. 

In this section we make use of xTensor [118], a free package for Mathematica [119], especially of its 
fast index canonicalizer based on the package xPerm [120]. 



40 



6 RES ULTS FOR HAMILTONIANS 



6.1 Post-Newtonian Expansion 

The post-Newtonian expansion of the ADM Hamiltonian has been well studied for nonspinning objects, 
for the second post-Newtonian level see [121], and up to and including the 3.5 post-Newtonian order see 
[34, 122-124]. From this expanded Hamiltonian the approximate equations of motion can be derived in 
a straightforward way. In this section we derive general formulas for the ADM Hamiltonian up to and 
including the formal second post-Newtonian order, which will then be applied to calculate spin corrections 
to the Hamiltonian in section 6.2. Another interesting application would be to obtain spin corrections to 
the post-Minkowskian Hamiltonian, see, e.g., [125] for the nonspinning case. 

Besides the ADM formalism, there are various other methods available for post-Newtonian calcula- 
tions. The equations of motion at the first post-Newtonian order are due to Einstein, Infeld, and Hoffmann 
[126], obtained with the help of a surface integral approach. This method got further developed and ap- 
plied up to and including the third post-Newtonian level, see, e.g., [127]. A further important method uses 
point-masses in harmonic gauge, which also succeeded to derive the third post-Newtonian order equations 
of motion; for a review see [5]. This method has advantages for flux and waveform calculations, which 
succeeded up to the third post-Newtonian order [128] (corresponding to the knowledge of the equations of 
motion at the 5.5 post-Newtonian level, which seem to be impossible to obtain directly). Another approach 
in the harmonic gauge is the direct integration of the relaxed Einstein equations, see, e.g., [129]. More re- 
cently also methods inspired by quantum field theory were developed, see, e.g, [110, 130]. An advantage of 
these methods is that some of the very sophisticated and systematic techniques for perturbative calculations 
used in high energy physics can be applied in a straightforward way. 



6.1.1 Review of the Formalism 

We now give a short summary of the calculation of the ADM Hamiltonian. First the field constraints 



167T^/7 
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n /matter ^ _ 1/ matter /£■ i \ 

= H , -— n/ij7r ', k =H i , (6.1) 



have to be solved within the ADM transverse traceless gauge, which for the metric leads to the decompo- 
sition 

7«= + ^+^ T ' (6 - 2) 

at least to linear order in spin. Such a solution can in general only be found in some approximation 
scheme and we consider the post-Newtonian one here. Having the decomposition (6.2) one can solve the 
Hamilton constraint for <j> (this will become obvious in the next section). Then we can calculate the ADM 
Hamiltonian 

Hadm = -72- / d 3 xA0, (6.3) 
16tt J 

which must be expressed in terms of the canonical variables. It is suitable to already express the source 
terms % mattel and %™ attel in terms of the canonical matter variables, which is done in section 6.2.1. Then 
no further redefinition of the matter variables is necessary. 

However, it seems to be simpler to perform the redefinition of the field momentum after solving the 
constraints. As the gauge condition at linear order in spin now reads 7r" = 0, or, with (4.18), (4.50), and 
BfjS kl = 0, 



ik jl mapakuSal 
1 7 7T— ( 

2np a {m a — n\ 

the decomposition (2.51) is not valid any more. But we can still use the general decomposition 



TT* = -lfor$>«f = -lforXXVV ., m f ~ J a , (6-4) 
^ *-f 2np a (m a - np a ) 



(6.5) 



6.1 Post-Newtonian Expansion 
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with 



r «TT 



7T • 



flltfl.-A- 1 ) 7T 



A A 



1 

2 J 



(6.6) 
(6.7) 



and the vector potential is still fr* = A - 7r* J j. This can be shown by inserting (6.6, 6.7) and (2.55) into 
(6.5), which then turns into an identity. The new part ■k %j can immediately be obtained using (6.4). After 
the constraints have been solved using this decomposition, we go over to the canonical field momentum 

Tt ijTT by 



r ijTT 



-ijTT 



(6.8) 



No redefinition of hJJ is needed at the linear order in spin. 
6.1.2 Expansion of the Constraints 

Now we expand the constraints according to the formal post-Newtonian counting rules introduced in sec- 
tion 4.2.3. Notice that only the field parts tfi, tt* 3 , and 7r* J are expanded, but not hJJ and w lJ . The latter 
are still dynamical variables in the ADM Hamiltonian and can be expanded only after their equations of 
motion were obtained and solved. For the Hamilton constraint we get 
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— W^-^j + (td) , 



(6.9) 



(6.10) 



(6.11) 



up to and including the formal second post-Newtonian order. These equations can be solved iteratively for 
<j> by applying an inverse Laplacian to them. The ADM Hamiltonian (6.3) results from an integration over 
the right-hand sides of these equations. It was used that 7r iJ = O (c~ 9 ) at linear order in spin. 

However, we also have to solve the momentum constraint as 7r u appears on the right-hand side of 
the Hamilton constraint. The expansion of the momentum constraint immediately follows from the exact 
formula 



71 ,3 = 



(6.12) 



with 



B 13 = 



l-(l + i0) 4 (7T« + 



1 



J T 3 k r 



r^ TT ) + V k (hH + hY,i - hjl k ) - -„V\ k hJ] , 



(6.13) 
(6.14) 



which is analogous to (4.21). Here we introduced the alternative vector potential 
1 



V 



(6.15) 
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for which it holds 

n ij = V\ j + V j , i -^S ij V k , k . (6.16) 
To the considered order we thus have 

*g> j - - 8 ^r > = -**>w - \ (^2)tf 3) ) . . (6.i7) 

With the help of tt % = A -1 tt 13 ■ , the expanded momentum constraint can be solved iteratively for 7p by 
applying an inverse Laplacian to it. tt' ij and V % then follow from (6.7) and (6.15). 



6.1.3 Formulas for Hamiltonians 

The first contribution to the ADM Hamiltonian (6.3) results from an integration over the first relation in 
(6.9) as 

H = J d 3 xH^ r . (6.18) 

Notice that %™" er is just the Newtonian mass density, so Ho is the constant energy belonging to the total 
Newtonian mass. Similarly, from the second relation in (6.9) we obtain the Newtonian Hamiltonian 



H N = / d 6 x 



nj matter ^ ± nj matter 
^(4) ~o^(2) tt (2) 



(6.19) 



( 2 ) results from (6.9) as 



■WttA^HT^ , (6.20) 



L (2) 

^•(2) ' "-(4) 



and agrees up to a factor with the Newtonian gravitational potential of the mass distribution H 1 }}^"- 'HJ nattcl 



is the Newtonian kinetic energy density. 

Next we proceed to the Hamiltonian at the first post-Newtonian order. However, we first apply the 
partial integration formulas 



- f,v. 1 
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to the right-hand side of (6.10), following from the formal solution 



i(-;3)) 2 =^)^ r +(td), 



(6.21) 
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matter 



(6.22) 



of the Hamilton constraint and from (6. 16, 6. 17). V? 3 n is determined by (6. 15) and 

4) = -8ttA- 1 ^. (6.23) 
Finally we get for the first post-Newtonian (PN) order Hamiltonian 

1 i .m'iHoi- 1 



•2 / matter " 1 nj matter ( " ±2 nj matter , jri nj matti 
W (6) ~ 4^(2)^(4) + 32^(2) H {2) + V (3) H (3)i 



(6.24) 



Notice that all terms in the Hamiltonian involve matter source terms % matter or Ti.™ and are thus integra- 
tions over delta distributions only. Further only the Newtonian potential (j>r 2 ) and the leading order vector 
potential need to be determined (lapse N and shift N l are not even needed at any higher order). This 
shows the efficiency of the ADM formalism in calculating the conservative post-Newtonian dynamics. 



6.1 Post-Newtonian Expansion 
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In the same way one can obtain a formula for the second post-Newtonian Hamiltonian 



i? 2 A P D N M = / d 3 x 
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where the partial integrations 
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= leTrtfoWST - -0 (2) ir" ) + (td) , (6.27) 



1 / l \ 2 

'(3)"(5) - iu " "(3)"-(6)i - 2^ (2) B)> 

were used. Notice that (f>^ and were eliminated from the Hamiltonian by these partial integrations. 
Therefore no solutions to the constraints besides (6.20), (6.22), and (6.23) have to be determined explicitly. 
The Hamiltonian i?2PN M nas tne additional label ADM as it still depends on the dynamical field variable 
hJJ. The elimination of hJJ from i?2PN M leads to the matter-only Hamiltonian i?2PN and is discussed in the 
next section. Further 7r ijTT first appears at the formal third post-Newtonian level. 

Notice that the obtained formulas are valid for quite general source expressions % mattel and %™ atter , not 
only to the ones linear in spin. 

6.1.4 Matter-Only Hamiltonian 

In the last section the ADM Hamiltonian Hadm was expanded as 

-Hadm = Ho + H N + iJiPN + #2pn M + • • • • (6.28) 

The conservative matter-only Hamiltonian results from plugging the solution for hJJ and 7r , -' TT into the 
action, Eq. (3.63) (and a subsequent elimination of emerging higher order time derivatives of the matter 
variables), see [123]. As tc 1 ^ is neglected at the considered order, the first term in (3.63) does not con- 
tribute here. Therefore -ff^™ turns into the matter-only Hamiltonian i?2PN by simply inserting the solution 

for^T into #ADM 

The field evolution can be obtained from the ADM Hamiltonian by 



dhij r ,TT tt i 1C rTTii ■ <5-Hadm 
-Qf - {fiij > adm j - lbird kl 5%km 



(6.29) 



{7r^ TT ,iJ ADM } = -167r<5™ J - 
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However, as the 7r yTT -contributions are of higher order here, we formally just have 

X IT ADM 

or explicitly, given that W^ ttcr has contributions linear in hJJ, 



5(J d 3 a;-H mattcr 



A//, 1 / = 2^ lfcI / (4)fei , with f m = 16tt v " ^ feTT (8) 1 ~ ^ (2 )^ (2)J . (6.31) 
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Using the formal solution h T J = 2SjJ kl A 1 f(4)kh a ^ contributions of hJJ to -£/2pn can be collected as 



If one is interested in the spin contribution of this integral only, one can obviously perform a partial inte- 
gration in (6.32) in such a way that only the spin part of hJJ is needed, see also [14]. This is desirable as 
the spin-dependent part of hJJ is much simpler than the spin-independent part. 

Though the discussion of hJJ was straightforward here, it is quite subtle to obtain the post-Newtonian 
expansion of (6.29) at higher orders. Indeed, it is not easy to correctly implement the boundary conditions 
into the solution of the first order equations (6.29). At higher orders (6.29) can be converted into a (second 
order) wave equation for hJJ, with source terms expanded according to the post-Newtonian counting rules. 
In [122, 123] this wave equation is then solved order by order using a near zone expansion of the retarded 
solution up to the 3.5 post-Newtonian order, corresponding to the boundary condition of no incoming 
gravitational waves; see also, e.g., [5] for other aspects like tails. Equation (6.31) is indeed the leading 
order near zone approximation of the wave equation for hJJ. The solution for hJJ at higher orders is 
responsible for the half post-Newtonian orders in the matter-only Hamiltonian, starting at the 2.5 post- 
Newtonian order. 

6.2 Spin Corrections to the Hamiltonian 

By now there are a lot of results regarding spin effects at the conservative orders in the post-Newtonian 
approximation. The main goal of this section is to derive the next-to-leading order spin effects within the 
developed formalism, which were tackled only recently. Even higher post-Newtonian orders linear in spin 
were derived recently in [44] for test spinning objects in the Kerr metric. Also Hamiltonians of cubic and 
higher order in spin were obtained for binary black holes [44, 116, 117]. The calculation of the leading 
order dissipative spin-orbit and spin(l)-spin(2) Hamiltonians was prepared in [15]. The corresponding 
equations of motion were already obtained [131]; see also the considerations in terms of orbital elements 



More work needs to be done for an application of the Hamiltonians derived in this section to gravi- 
tational wave astronomy. In particular the spin contributions to the next-to-leading order radiation field 
are only known for the spin-orbit case [102], but not yet for the spin(l)-spin(2) and spin(l)-spin(l) cases 
(for the latter case the stress-energy tensor derived in section 5 is needed). Further, it would be useful to 
find a parametrization of the orbits by solving the equations of motion, i.e., extending the solutions from 
[133, 134] at least to some of the new Hamiltonians. Finally, one should consider to incorporate the new 
Hamiltonians into the very successful effective one-body approach [9], which already succeeded for the 
leading order spin Hamiltonians [10] as well as for the next-to-leading order spin-orbit Hamiltonian [11]. 

Formulas and regularization procedures for the integrals that need to be solved in this section are given 
in, e.g., [122, 123, 135]. Some parts of the calculations were also checked using Riesz kernels in arbitrary 
dimension, see, e.g., [35]. 

6.2.1 Field Constraints in Canonical Variables 

Before starting the calculation of the Hamiltonians, it is suitable to express the source terms of the con- 
straints % matter and %™ attel " in terms of the canonical matter variables. Then the formulas provided in 
section 6.1.3 automatically give the Hamiltonian (the redefinition of 7r ijTT is not necessary here). Applying 
the variable redefinitions from section 3.2.2 to (3.39, 3.40) leads to 
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in [132]. 
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(6.34) 
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where 



S a = 7' 



<->mfc +7 7 — ~i 

np a (m a - np a ) 



(6.35) 



and A M given by (3.44). These source expressions are valid in general, also within the test-spin Hamilto- 
nian (3.47). In the spatial symmetric gauge it holds A 1 ^ e^ n ) k e^ m \^ — Tr^ lr fkl,n> where (3.55) and (4.49) 
were used and 71™ is given by (4.50). Notice that the variable redefinitions from the action approach lead- 
ing to these expressions have been checked up to and including the formal 3.5 post-Newtonian order by 
the symmetry generator approach [15]. This includes the formal third post-Newtonian or next-to-next-to- 
leading order linear in spin, which for maximal spin is at the 3.5 post-Newtonian order in the spin-orbit case 
and at the fourth post-Newtonian order in the spin(l)-spin(2) case. It was shown in [15] as a further check 
to the same level of approximation that the wave equation for hJJ following from the ADM Hamiltonian 
agrees with the Einstein equations, which again verifies that the used variables are canonical. 
The expansion of % matter sufficient for the formal second post-Newtonian Hamiltonian reads 
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where p a = (p a i)- The expansion of the source %™ atter is given by (6.34) and 
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Notice that the triad terms in (6.34) do not contribute at the considered order. 

The expansion of the static source terms needed at the spin(l)-spin(l) order follow from (5.34) as 
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and it holds 
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Here S a = {S a u\) and S a u\ — ^£ijkSa(j)(k)- No further contributions to %™ atter arise in the spin(l)- 
spin(l) case. 
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6.2.2 Leading Order 

The leading order spin effects are at the formal first post-Newtonian order and their Hamiltonian can be 
obtained from (6.24), which of course gives the first post-Newtonian Hamiltonian in the nonspinning case. 
The needed solutions of the constraints read 
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where f a — |x — z a | and z a — (z a ). The leading order (LO) spin-orbit (SO) Hamiltonian follows as 
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where f ao = |z a — z&| and n a & = (z a — Zf,)/f ao . This Hamiltonian is at the 1.5 post-Newtonian order for 
maximal spins. Further, the leading order spin(a)-spin(6), or S a S , Hamiltonian results as 



H s?s b = E E 

a b^a 



1 

2f? 



3(S Q • n ab )(S b ■ h ab ) - (S • S 6 ) 



(6.47) 



For maximal spins this Hamiltonian is at the second post-Newtonian level. Finally, the leading order 
spin(a)-spin(a), or S Q , Hamiltonian is given by 
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which is also at the second post-Newtonian order for maximal spins. All Hamiltonians in this section are 
valid for arbitrary many spinning objects. The Poisson brackets are the standard canonical ones, i.e., 



(6.49) 



zero otherwise. 

The leading order spin effects derived here are well-known for black holes (Cq — 1), see, e.g., [136, 
137]. For the leading order Cg-dependence see [105, 137]. 



6.2.3 Next-to-Leading Order 



Now we proceed to the formal second post-Newtonian Hamiltonian (6.25), which includes the next-to- 
leading order spin effects. There we also need the functions 
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Notice that there are CQ a -contributions in (6.50). We restrict to two spinning objects in this section. The 
results provided here complete the knowledge of spin corrections to the Hamiltonian up to and including 
the third post-Newtonian order for maximal spins. 
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Next-to-Leading Order Spin-Orbit Following the method developed here, the next-to-leading order 
(NLO) spin-orbit Hamiltonian results as [14] 
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(6.52) 



where (1 f-> 2) indicates an exchange of particle labels, and is identical to the one derived earlier in 
[84]. The next-to-leading order spin-orbit case was first tackled on the level of the equations of motion 
in [138] and was later rederived and improved in [139] (both in the harmonic gauge). Within the ADM 
canonical formalism the Hamiltonian H^q° corresponding to these equations of motion was obtained in 
[84] from the spin equation of motion (2.28). The linear-in-G part of H^q was also derived in [116] 
from corresponding source terms of the constraints, similar to the approach used here (however, in [116] 
the source terms were obtained from the approximate Kerr metric in the ADM transverse traceless gauge). 
Very recently derivations within the effective field theory approach also succeeded [140]. 



Next-to-Leading Order Spin(l)-Spin(2) The spin(l)-spin(2), or S1S2, Hamiltonian reads [20] 
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and was confirmed by [114, 141]. Notice that no agreement with the result in [142] could be found, see 
[20]. Indeed, the result in [142] turned out to be incomplete [20, 1 14]. 



Next-to-Leading Order Spin(l)-Spin(l) A nonreduced potential (i.e., with the spin supplementary con- 
dition not eliminated on the level of the potential) for the next-to-leading order spin(l)-spin(l), or S^, 
dynamics is given in [107, 143]. Within the method described in section 5.3.2 an equivalent Hamiltonian 
i/NLO w jjj k e (j er i vec j here. This Hamiltonian was first given only for the black hole case (Cqi = 1) in 

[17] and then generalized to arbitrary Cqi later [19]. However, the comparison with [107, 143] was quite 
cumbersome. First agreement with [107] could not even be found in the spin precession equation [17], 
however, this finally succeeded after identifying a sign typo in [107], see [18] (all for the case Cqi = 1). 
After a further correction [143] full agreement was finally found in [19], now also for arbitrary Cqi. For 
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this comparison the potential from [107, 143] was first transformed into a fully reduced Hamiltonian in [19] 
by a Legendre transformation and an elimination of the spin supplementary condition using Dirac brackets 
(2.75). Then a canonical transformation leading to our result in [19] was searched for and found. 
The result for general compact objects (including neutron stars) is [19] 
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The corresponding spin(2)-spin(2) Hamiltonian simply results from an exchange of particle labels. 

According to section 5.3.2 the linear-in-G part was derived with the help of the Poincare algebra method 
from [117], while the G 2 part (the last two lines) results from the source expressions (6.40-6.42) derived 
in the present thesis. 

Notice that for black holes (Cqi = 1) this Hamiltonian was already found in [17], for the first time 
including the correct center of mass motion. Further, the earlier result for the general case in [107, 143] 
is not a fully reduced Hamiltonian. The Hamiltonian presented here is on a higher level of sophistication 
with advantages for applications, e.g., the spin vectors appearing in our Hamiltonian have a constant length 
and it is easier to obtain all equations of motion in terms of these "good" spin variables. 



6.2.4 Center of Mass and Poincare Algebra 

The post-Newtonian expansion of the center of mass vector 
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can be obtained from the expanded Hamilton constraint (6.9, 6.10). To the formal second post-Newtonian 
order this leads to 
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For the formula for G, PN partial integrations were applied, similar as in section 6.1.3. Notice that tt^ k is 
spin-independent. For results in the nonspinning case see [124]. 
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The contributions to the center of mass vector corresponding to the leading order spin Hamiltonians 
follow from Gipn = {G\ PN ) as 
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From G2PN = (G\p N ) the next-to-leading order parts result as 
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Notice that (6.59) and (6.60) are valid for arbitrary many spinning objects, while (6.61) holds for two 
objects only. For two objects GgQ° was already found in [84]. Further, G^ simply results from an 
exchange of particle labels in (6.61) 

Now one can check whether the Poincare algebra (2.10, 2. 1 1) is fulfilled, which is indeed the case (the 
Hamiltonian plays of course the role of the energy E). At the spin-orbit level this was already shown in 
[84]. At the spin(l)-spin(l) level this holds by construction, as most terms of the Hamiltonian H^t° were 
obtained from the Poincare algebra via an ansatz in [1 17]. However, the fulfillment of the Poincare algebra 
provides a thorough check of H$q and 



7 Conclusions and Outlook 

The first main goal of this thesis, the extension of the ADM canonical formalism from nonspinning to 
spinning objects, succeeded to linear order in spin via an action approach. The result was verified by 
an independent order-by-order derivation. Even the extension to higher orders in spin is well understood 
now, but somewhat more complicated and requires further approximations, like the post-Newtonian one. 
The second main goal of this thesis, the calculation of conservative Hamiltonians for inspiralling binaries 
relevant for gravitational wave astronomy, was then straightforward. The effort of first deriving the canon- 
ical formalism was payed off by its efficiency in the calculation of these Hamiltonians. New results are 
the next-to-leading order spin(l)-spin(2) and spin(l)-spin(l) Hamiltonians, and the spin-orbit Hamiltonian 
derived earlier by Damour, Jaranowski, and Schafer was confirmed. All Hamiltonians through the third 
post-Newtonian order for maximal spin are known. 

The next most interesting Hamiltonian which could be calculated is the conservative next-to-next-to- 
leading order spin-orbit one, which is at the 3.5 post-Newtonian level for maximally rotating objects. No- 
tice that the verification of the canonical formalism given in this thesis via the order-by-order construction 
already covered this case. Leading order dissipative Hamiltonians are also envisaged and its calculation 
was already prepared in [15]. For maximal spins these Hamiltonians are even at the fourth post-Newtonian 
order in the spin-orbit case and at the 4.5 post-Newtonian order in the spin(l)-spin(2) case. The exten- 
sion of a recent result within the post-Minkowskian approximation [125] to spinning objects would also 
be desirable, as it could be applied to the gravitational scattering of spinning bodies moving at relativistic 
speed. 
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Further, more work needs to be done for an application of the new Hamiltonians presented in this thesis 
to gravitational wave astronomy. In particular the spin contributions to the next-to-leading order radiation 
field are only known for the spin-orbit case [102]. This result should be extended to the spin(l)-spin(2) 
case, as well as to the spin(l)-spin(l) case. For the latter the spin(l)-spin(l) contributions to the stress- 
energy tensor given in this thesis are crucial. Also an implementation of the new results given here into the 
very successful effective one-body approach would be appealing. 

Another though rather mathematical development for the future would be to consider the full constraint 
algebra, gravitational field and supplementary conditions, at different stages of gauge fixing, as well as a 
treatment using Dirac brackets; see also [64] for the case of Dirac fields. 
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Symbols 






A* 


defined by (3.44) 


gij 


defined as e ij = \ {e 1 ^ ~ e j ^) 


D kl 


defined by (3.54) 


gij 


defined as e l J = e iJ (z^) 


c 


speed of light, usually c = 1 here 




3-dimensional Levi-Civita symbol 


Cq 


mass-quadrupole parameter, see (5.11) 


E 


energy of the system, see (4.19) 


S 


defined as 5 — 5(x l — z l ) 




electric part of cj£ ap , see (5. 19, 5. 18) 


5 


defined as 5 — 5(x l — z l ) 


F 


time like vector in conditions (3.12, 3.10) 


<5(4) 


defined as (5( 4 ) = 8{x^ — z^) 


<t> 


trace part of the induced metric, see (6.2) 


Sij 


Kronecker delta, (Sij) = diag(l, 1, 1) 


9 


defined as g = de^g^) 


SjTTkl 

ij 


transverse traceless projector, see (2.55) 


SV 


4-dimensional metric 


A 


Laplace operator, A = didi 


7 


defined as 7 = det(7y-) 


A- 1 


inverse of A for usual boundary conditions 


lixv 


projector (2.38), contains induced metric 




tetrad field, g^ v = e/ M e J v 


G 


gravitational constant, usually G — 1 here 




triad in the symmetric gauge, see (3.53) 


G l 


center of mass vector, see (4.19) 
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Tkij 3-dim. Christoffel symbol of first kind 

ri 4 /L 4-dim. Christoffel symbol of first kind 

hJJ transverse traceless part of 7^, see (6.2) 

H general symbol for a Hamiltonian 

^adm ADM Hamiltonian, see (3.61) 

H Hamilton constraint, see (6.1) 

Hi momentum constraint, see (6.1) 

H md field part of H, see (2.44) 

Hf ld field part of Hu see (2.44) 

^matter matter part of H, see section 6.2. 1 

^matter ma ^ eT f see section 6.2.1 
^matter defined by (4 2 2) 

I moment of inertia of a spherical top 

J Miy total angular momentum, see (2.12, 4.9) 

jvap Dixon's quadrupole moment 

Kij extrinsic curvature, see (2.39) 

A Lagrange multiplier 

Lm matter Lagrangian 

Cg held Lagrangian density (2.36, 2.40, 2.43) 

Cm matter Lagrangian density 

A AI Lorentz matrix, see section 3.1.3 

AM W canonical rotation matrix defined by (3.43) 

to constant mass-like parameter, see (3.28) 

too constant mass-like parameter, see (3.26) 

to p dynamical mass, = — 

M mass of the system, M 2 = -P^P^ 

normal vector for (3+l)-split, see (2.37) 

h ab defined by n ah = (z a - z b )/f ab 

np defined by np — n^p^, see (3.32) or (5.30) 

np defined by np — — J m 2 + j^Pipj 

nSi defined by nS t = n^S^, see (3.33, 3.42) 

N lapse function, see (2.37) 

N* shift vector, see (2.37) 

uj^ ij Ricci rotation coefficients 

fl fiV angular velocity tensor, see (3.15) 

angular velocity fiWO') = A [k] ^kW^ 

p^ linear momentum, see (2.31, 5.4) 

Pi canonical momentum conjugate to z l 

P^ total linear momentum, Pq = —E, (4.9) 

P M „ the projector P^ v = - j^pf^f" 

n ij defined by (2.41) 



Ti-y'TT transverse traceless part of 7r u , (6.5, 6.6) 

tt 1 J spin correction to 7r 4:,TT , see (4.18) 

tt 1 vector potential for 7T y , see (6.7) 

ir 4 - 7 vector potential part of 7r u , see (6.5, 6.7) 

tt^ trace part of ir^ , see (6.5, 6.6) 

■fr^ canonical field momentum, see (3.58) 

Tj-y'TT transverse traceless part of tt^, see (3.60) 

7T 2 vector potential for ir 1 - 7 , see (3.60) 

tt 1 ^ longitudinal part of 7r y , see (3.60) 

Q^ v mass quadrupole part of J^ ua P , see (5.6) 

Qij defined by (5.35) 

f a defined by f a = |x - z a | 

f a6 defined by f a6 = \z a - z 6 | 

i? 3-dimensional Ricci scalar, i? = 7* 3 'i?y 

i?( 4 ) 4-dimensional Ricci scalar, i? = gP"Rj$ 

Rij 3-dim. Ricci tensor, i?^ = J Rikjl 

R$ 4-dim. Ricci tensor, = ff^-R^ 

3- dim. Riemann tensor, sign as in (2.29) 
R^l a p 4-dim. Riemann tensor, see (2.29) 

4 j defined by (4.30), see also (6.35) 

S spin length, 2S 2 = S^S^ 

spin tensor, usually restricted to (2.32) 

S^ v canonical spin tensor, see (3.49) 

Sa) canonical spin vector, Sa\ = \^ijkS(j)(k) 

S canonical spin vector, S = {Su\) 

t time coordinate, a; = t or z° = t 

t worldline parameter 

t> iV — multipole moments, see (2.26) 

S9 IJ variation for A AI , 56 IJ = A j4 / M AJ 

T^ v stress-energy tensor 

(td) denotes a total divergence 

4- velocity, = ^ 

V 1 vector potential for ir 4 - 7 , see (6.15, 6.16) 

V 1 vector potential for , see (4.24, 4.25) 

W full action, W — Wq + W M 

Wq Einstein-Hilbert action, see (2.36) 

Wm matter part of the action W 

x M spacetime coordinates, x° = t 

z^ worldline function, z° = t 

z 1 canonical position variable 

z\ possible correction to z l , see (4.34) 
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